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S YNOPSIS 

MOSSBAUEP STUDIES OE THE HYPERTINE INTEHACTIOIS 
IN SOHB hep METALS AND DISORDERED BINARY ALLOYS 

By 

N. VASUDBVAN NAIR 
DEPARTMENT OP PHYSICS 
INDIAN INSTITUTE OP TECHNOLOGY 
KANPUR 208016, INDIA. 

MUssbauer effect, v<hich is the recoilless absorption 
or emission of gamma rays by nuclei, has provided a powerful 
technique for studying the interactions of the atomic nucleus 
with its environment. The interaction of the electric quadrupole 
moment of the nucleus and the electric field gradients (EPG) 
present at the nuclear site is called the electric cuadrupole 
interaction, and that of the nuclear magnetic dipole moment 
with the magnetic field at the nucleus is known as the magnetic 
hyperfine interaction. 

It is found that in many of the non-cubic metals the 

electronic contribution to the EPG is about 2-4 times the 

lattice contribution. It is also interesting to note that 

with all the complexities, the dependence of EPG on temperature 

■5/2 

obeys a simple T"^'' relationship. More experimental data on 


the non-cubic metals will be valuable to find the exact syste- 
matics and to improve the theoretical understanding of the 



origin of EPG-. The present thesis, in part, deals with the 
measurement of EPG- and its dependence on temperature in some non- 
cubic systems. These results are analysed on the basis of some 
of the available theoretical models. 

Mbssbauer spectroscopy has been found verjr useful in the 
investigations of order-disorder transitions in alloys and its 
effect on hyperfine magnetic field and quadrupole interaction. 

The broad distribution of magnetic hyperfine fields in the 
disordered alloys is notably reduced by atomic ordering. Metal 
atoms coagulate some times into the non-crystalline form which 
is in an extraordinarily disordered state called the amorphous 
state. The amorphous rare-earth based alloys have raised conside- 
rable interest because of the possibility of its use in computer 
memories and thermomagnetic recording. The magnetic properties of 
these alloys are quite varied and sensitive to the variation of 
temperature and other parameters. The study of an amorphous 
system leads to a critical re-examination of the role of lattice 
in crystalline solids in the structural origin of hyperfine 
fields. The Mbssbauer and magnetic studies of two such 
disordered binary alloys are reported in this thesis. 

The methods of analysis (BEG- distribution and recoilless 
fraction) developed by us were first checked on certain standard 
systems, e.g, Ee(SO^ ).7H20, A thorough systematic study of these 
old systems yielded some outstanding new results such as identific 
tion of the induced Stark effect, which is presented briefly in 
this thesis, 

A general brief introduction to Mbssbauer spectroscopy 
is given in Chapter 1, Mbssbauer parameters like the recoilless 





fraction, isomer shift, electric quadrupole interaction, 
magnetic dipole interaction are introduced. Other topics like 
(a) Relaxation effects (b) G-oldanskii-Karyagin effect (c) 
Sternheimer antishielding effect and (d) Asymnietric doublets 
in Mossbauer spectra are discussed briefly in this chapter. 

Chapter 2 gives the experimental details including the 
authors contribution in this field, e.g. designing and develop- 
ment of a ‘miniature quartz vacuum furnace (Induction heating/ 
resistance heating) for H'dssbauer studies'. The method of 
preparation of the samples and the optimisation of the thickness 
for best signal to noise ratio are discussed. The computer 
programmes used in the analysis of data are also briefly 
explained. 

Chapter 3 presents some new studies done on some well- 
known ferrous compounds of iron. The effect on ERG- by an 
externally applied electric field to these samples is described 
at the first part of this chapter. The later port gives the 
study of the ERG distribution in these compounds by analysing 
the Mossbauer data by a computer programme developed by the 
author using a model independent method. The recoilless 
fraction and the Debye- temperature of some of the ferrous 
compounds are determined. 

In Chapter 4, the present status of the understanding 
on the origin of EPG's in pure metals and alloys is briefly 
reviewed. The systematic trends observed in the experimental 
findings are discussed. The theoretical models suggested 
for finding the EEG in hexagonal close packed (hep) metals 
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are also reviewed? especially the charge screening approach of 

Mshiyama et al., and the change shift model suggested hy 

57 57 

Bodenstedt and Perchied, The measurements on PeTi , PeSe 

57 

and FejSTd. systems at various temperatures are reported and 
discussed in detail. 

The quadrupole interaction frequency (QIP)? at Pe as 

a dilute impurity in Tij,as an absorber in the M'dssbauer 

experiments is determined for the first time. The measurements 

are done from 80 £-586 £, The dependence of BPG on temperature 

obeys the T^^^ law as for most of the other systems. The lattice 

contribution to PPG- is calculated using the Das and Pomerantz 

relation and the electronic contribution is separa,ted from the 

experimentally observed EPG. The present measuiement is in 

agreement with the correlation observed between the electronic 

57 

and ionic contribution. The PPG' a at Pe in Ti at different 

temperatures are calculated using the charge shift model aiid 

it is found that there is a general agreement between the 

theoretical and experimental values. 

57 

The SPG at Pe probe nuclei in the semiconductor 

selenium is measured in the temperature range 82 £ - 298 £, 

57 

The present data on "PeSe seem to qualitatively support the 
correlation between the ionic and electronic contribution to 
EPG if we assume that the experimental and ionic EPG's are of 
opposite signs. Otherwise the electronic contribution to the 
total PPG becomes negligible compared to the ionic contribution. 
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5 7 

The recoilless fraction of FeSe absorber is measured 

using the method of Tr coster et al. with some modifications. 

57 

The EFG- at Fe in the light rare-earth metal neodymium 
is measured for the first time in the tempera'ture range 82 K - 
295 K. It is found that this system violates the systematics 
proposed by Ernst et al. regarding the ratio of the electronic 
to ionic contribution to EFG. It is also found that the 
electronic contribution to the EFG is the main source of EFG 
in this system. 

In Chapter 5 the Mbssbauer and magnetisafcion studies 
done on the Fe-31.2 /, at Pd alloy and the Mq Fe^ 
amorphous rare-earth alloy are reported. The pha,se transition 
from fee to fet at around 173 K is observed in the Fe-Pd 
system. The systematic behaviour of the quadrujpole interaction 

3/2 

with temperature in the fet phase obeys the T ' law. The 
Mbssbauer data are analysed by the hyperfine field distribu- 
tion method developed by Window and the width of the field 
distribution is found to increase with temperature, which 
reveals that the range of exchange interaction is short in 
the Fe-Pd system. The universal correlation that the total 
EFG is proportional to the lattice contribution is obeyed by the 
tetragonal phase of Fe-Pd .alloy. The recoilless fraction of 
this alloy at various temperatures are determined by the 
Trooster et al. method. The paramagnetic Curie temperature 
is determined by using the PARC vibrating sample magnetometer 
and also from the variation of 


with temperature and the 
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agreement between these t'wo values is quite good in this case. 

The magnetisation studies done on the amorphous system 

NIq Fe^ is presented in detail in this chapter. The 

Curie temperatiire is found out from the Mbssbauer studies 

and from the Arrott plots. The average hyperfine field varies 

from 317 kOe to 109 icOe in the temperature range 90 X to 320 E, 

The hyperfine field distribution is also found in this 

temperature range. In the paramagnetic phase this system shows 

quadrupole splitting with well defined peaks. The quadrupole 

■5/2 

splitting changes with temperature according 'co tbe T ' law. 
The distribution of EFG in the paramagnetic phase is also 
determined in the Mq ISq system. 

It is hoped that the present work has contributed to 
our knowledge on the systematics of hyperfine interactions 
in metals 5 alloys and compounds and has helued bring us closer 
to the understanding of the origin of the hyperfine fields. 



GItAPTSE 1 


EUNDAKEETAIS of H0S33AOER 3PEaTEOSGO?Y 
1,1 Introduction 

The discovery by R.L. Mossbauer [1] of the recoilless 
absorption or emission of gamma rays by nuclei has opened up 
new vistas in the realm of riyperfine interactions. The 
unlimited potentialities of this discovery is reflected in its 
rapid development into a scientific discipline and the 
impressive position it has occupied in the varied fields of 
physics, chemistry, metalliurgy, biochemistry, earth sciences 
to name a few. It is exciting to realise that by rxo'vi a single 
aspect of M'dssbauer effect like the isomor shift could fill 
the whole volume of a book [2], 

The field of hyper fine interactions with special 
reference to M'dssbauer spectroscopy is treated in many reviews 
and monographs in details Therefore only a broad outline 

will be given of the nuclear and solid state properties that 
influence M’dssbauer spectra in thns chapter, 

1.2 The M'dssbauer Effect 

When a nucleus decays from an excited state of energy 

1q with respect to the ground state of energy it will 

emit a photon of energy 


To conserve the momentum a 



( 1 . 1 ) 


3C0il energy Bj^ is imparted to the emitting nucleus 

,2 

Ji 

Er 




E 

w — w I 


2Mc 


where M represents the mass of the nucleus and c the velocitj' 
of light. Conservation of energy yields for the photon energy 

— 0X 

B = Eq-E^ similarlj/ to excite such a nucleus from B^ to ii > 
a photon of energy B = Eq + E^ will he required, it is clear 
therefore that resonant absorption of the emitted photon can 
only take place ^ when Bj^ is comparable to oi" smaller than r , 
the line width (full width at half maximum or in short BWHI4) 
of the excited state B'^'^. This line width is connected to Ihe 
mean life m of E®^ hy the Heisenberg uncertainty relation 


T 




Jl 

2 % 


where h is the Plank’s constant. In the case of gamma ray 
transitions, the recoil energy is much larger {’^ 10 times), 
when compared with typical nuclour line viidrhs thus making 
resonance absorption almost impossible. Ilow.jvcr . 'when the 
nuclei are embedded in a crystal, a fraction *f' oi the gamma 
rays are emitted without any energy loss duo to recoil, the 
same being true for absorption also. The probability for this 
is given by 

f = j/<p*(x) fpjj(x) exp (-i k.x) dx i ^ (1.2) 

where ^j^Cx) represents the vibrational state of the lattice, 
which remains the same in the case of recoilless 'emission, k 
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ivS the wave vector of the ga-^ma ray and x i? the position coordi 
nate of the centre of nasc of the omitting nucleus. 

A general expression for the probability of omission or 
absorption w'ithout recoil, that is in which no phonos are 
excited is given by 


f = exp <x‘^> ] 


exp [- 


v,2, 2 

h k 
2 m 


r Goth (T^rr™pr) dfi) J 

1 h 


<x“> is the 'mean square displacement of the nucleus from its 
equilibx’ium position, and g( “ ) is the normalised density of 
phonon states of angular froquency u , In the Debye model of 
solids, the frequency distribution function 


g (u) = const art , 


2 

oa 


f 0 r 0 ^ ^ 1 . 


max 


= 0 


for t»-' ^ - 

' — i. 


Using this model the oxpr oscion foi* 


I f 1 v- 


reduces to 


f = exp { 


O.IJ-T J 

ml 


! rri 9 . 

key h ! 




0 


X dx 
e’--l‘ 


1 } 


(1.3) 


It is seen that ’ f ’ con bi, non zero rindcr the- conditions 
mentioned above. The energy distribution of rocoilfree emitted 
radiation can be shown to be 

,2 


1(B) 


0 "^s ^ 


I 


r^/4 V 


(E-E^)‘^ Ar/4 


(1.4) 


while the probability for resonant absorption can similarly be 
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d escribed by 


c;(E) == 



( 1 . 5 ) 


whei’e and are the recoilfree fractions in source and 
aDsorber respectively;. is intensity of emitted radiation 

from the source and is the effective photon cross section 
of the Mbssbauer nucleus in the absorber. This tyne of recoil- 
less resonant absorption of nuclear gamma radis^tions is knot;n 
as Mbssbauer effect (ME) after its discoverer. 

Mossbauer effect has provided us with an extremely 
precise method for measuring the small en'srgy shifts experie- 
nced by the nuclei such as the hyperfine structure of tiie 
nuclear energy levels. 

In Mbssbauer experiments either the source or the 
absorber is usually subjected to a velocity programiae to give 
an appropriate Doppler energy shift to the Y--ays so that one 
may obtain the resonanco intensity as a, function of energy. 

1 . 5 flyporfino Interactions 


When a nucleus is .jmbedded in a solid, the Hamiltonian 
for interaction between the nucleus and the surrounding charges and 
currents may be written as 

A- + Ml + ^2 + 

where refers to electric mono pole interaction between the 
nucleus and electrons, the magnetic dipole interactions and. 
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Ep to electric quacrrpole interactions. Higher terris are 
nsually small. 

The electric monopoly intera.ction between the electrons 
and nucleus a hi its the eae:rrj j_evols of the gror.nd and the 
excited states of ■’■•h''- nucleus generally by cli7rf extent amounts. 
Depending on the chexiical foom, this change- ma;-- be different 
for source and absorber whren will cause a shift in the 
resonance energy. This o>j±j;fc is known as isoraei’ shift or 
chemical!- shift. 

The electric quadrupole and magnetic interactions 
result in splitting of the energy levels and consequently 
give a multiple resonance absorption pattern in the Mbssb-auer 
spectra. These interactions will be discussed briefly. 

1.5.1 Isomer shift 

The isomer sl'i"'‘t 5 j s the result of "^ho Coulomb inter- 
action of the nuclear charge uistribution because of the 
finite size of the nucleus th the electron charge density at _ 
the nucleus (s electrons). Because of the different environment 
of the source and absorber nuclei this interaction energy may 
be also different. So in a resonance experiment the centre of 
gravity of the spectrum will be shifted by this difference of 
interaction energies. 

This shift can be expressed in non-relativistic approxi- 
mation, assuming a spherical uniformely charged nucleus, 'as 
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6 


2 % 

T 




<r^> 

e 




( 1 . 6 ) 


vriiere the suffixes A and S rr::i7?er to absorber and source 

respectively. cp„(0) is the s electron density at the site 

2 2 

of nucleus and <r > and <r‘Z > are the mean squared nuclear 
charge radii in the excited and ground states respectively. 
(Pig. 1.1 ). The vreal th of information that could be obtained 
from isomer shift measurements are described in the recent '^rork 
on isomer shifts edited by l.ii. Shenoy et al. [2]. 

The ratio of isomor shifts of two gan:ma- transit ions A 
and 3 in two isotopes of the same element is 


r 2 

I <r > 

^ s . -A 


<rrr> 1 


in 


r 2 
[ <rg> 


5 ^ 


B 


<r > I 

S ^ 




which is independent of the electron density. That is, this 
ratio of isomer shifts gives the ratios of the changes in the 
mean squared charge radii for different transitions in the 
same element. 

The isomer shift data are commonly used to study the 
chemical structure, bonding of valence states etc. in compounds. 
The measured isomer shift values of the same element in different 
systems are related to changes in the electron densities and are 
therefore, used to find the structure and phase of the systems. 




^ . SOURCE (S) 

: ^ ; - . , , ,6> ;2K/s:^2e 0»| ;-t ;i(jsC03 J I ) 4 K r| 

i I i '•■■ I ; i n ; r f ■; i ;,; Jj; ; 

EJg.M Effect of; .monopple' tnteroct'ion; joni. the ^energy jl«vg|<, of 
absorbeH.'" 


ABSORBER (A), 



1.3.2 Electric quadrupole interaction 


The electric qnadrujrole interaction iaiailtonian is 
given by 


^0 = 




( 1 . 


where 


3^v 



3Xj gx^. 




= 

I 

e 



= / d^x p(x) 

tJ 



The Hamiltonian can be written in the terms of the 
irreducible comijonents of Q and v 

-. 2 -, 


Qo = 

eQ 

[7I3 

■?l{^2i-i") 

# = 

eQ 

■f 6 

' 

21 ( 21-1 ) 


At2 

©Q 

f 6 

Q2 = 

21(21-1) 

-r- 


(I,)' 


V = 
o 




2 ^zz 


+ J:^ (Y + iy ) 
^ y-g ^ xz - yz^ 


y 


+2 


(y - y + 2i V ) 

2^5 ^ XX Vy - xy^ 


where 
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Q = 


1 


<II 


i > 




3 "9' 


On transformation of oh© co tensors to the set of principal, axes 
of electric field gradient (BPS) tensor, the irreducible 
components of become 


V 


^ 'Iz = i SI 


'Vi= ° 


^2 = 


aye ^hx"^yy^ " JyS 


peq 


The quantities p and q are defined by the above equations. 
The Hamiltonian then becomes 



I 


q: 


TO 



e^Qa 


41 (.21-1) 


- I‘ 






If)] 


( 2 . 8 ) 


The separation of nuclear energy levels is rsually much 
larger than the magnitude of 'Kq* So, to a good approximation, 
the perturbation calculations require only the matrix elements 
of the form <Im| > which are diagonal in all quantum 

numbers except m. 

I 

Por axially symmetric EFG (p = 0) the matrix of this 
Hamiltonian is diagonal in jlm > representation. The energy 
shift in the state ll+m> for I > 1/2 is 



LTilOJ t>|C\J 


10 


e^Qq 

41(21-1) 


[3J 


0 




1 


e ^Qg 

“a 


which for I = n = + | is and I = m = + ^ is 

There is no shift in energ;/ of the states having I = 

Thus the quadrupole split spectrum is a doublet (Fig. 1.2), 
The situation may be complicated if the asymmetry parameter p 
is non zero. The matrix for is not diagonal in ! Im > 

representation and the off dia.gonal terms are 


<Im’ I 1 Im> = [(I+m-l)(l+m)(l+m+l)(l+m+2)]^.6j^, 

e^ O g 

where A = 4'j 1[ "j ' lp" " ) ' • This causes mixing between states having 

Am = 2. The secular eouation for some of the half integral 
values of I are given below [8]. 


Seciilar equation 
(B/A)^ - - 9 ^ 0 


(E/2A)^-(3n^+2l)(B/2A)^ - 20(l-p^) = 0 

(E/3A)'^-42(l+p^/3) (B/3A)^ - 64(l-p^) + 105(l+p^/3)^ = 0 

(E/6a)^ - 33(1+p^/3)CB/6a)^-44(1-p^)(E/6A)^ 

+ 44 { 1+p /3) (E/6A) + 144 (1+ p^/3)(l— p^) = 


For I = 'g the exact solution is 

1 5^2 = 3A (1 + n2/3)^ 

B ^^2 = U + 


(2.9) 


0 
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For other values of nuclear spin, the secul.ar equation is solved 
numerically. 

The quadrupolfc interaction does not remove the (+ m) 
degeneracy. 

1.3.3 Sternheimer antishi elding factors 


The strength of the quadrupole interaction at the 
nuclear site will he considerably affected due to the presence 
of the atomic electrons. The effective EPG a.t the nuclear 
site is usually written as 


-y-nucl 

^22 


(1 - yP V 


ext 

Z2 


+ (l-.Il) V 


1 00 oil 


ext 

where V__ is the SPG at the nuclear site due to the ions of 

the lattice in the absence of atomic electrons. An additional 

field (-Y . ^e^t^ results at the nuclear site due to the 

CC 22 

polarisation of tne el^ctxcn cloud. The v:-l---;3 of (I-Y--) 
the medium and heavy ions ore in the range 10-80 showing that 
the f: eld at the nucleus is considera,bly enhanced and therefore 
(1 -Yq^) is Called the antisiii elding factor. If the field sources 
are not entirely external, as in the case of rare-earth metals 
or in the case of ions having non-spherical outer shells, 
significant local contributions to tho EPG are present. The 
factor R corresponds to the atomic shielding of the EPG at the 
nucleus due to the unfilled orbital The values of R 

range fipm - 0.2 to -t 0,2, The above picture is not valid for 
metals and alloys because the IPG contribution at the nuclear site 
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cannot be pictured '‘.'holly at? due to local and non-local sources 
and it is well established that very significant contributions 
como from the conduction electrons. 

1,5.4 Magnetic dipole interaction 

A nucXeus witii spin > 0 has a magnetic dipole moment 
p, which can interact with i magnetic field ii at the nucle'us, 
resulting in the splitting of energy levels (i-Juclear Zeeman effect). 
This interaction Hajmiltonian 

-p.H = - i.H (2.9) 

G "ii 

whore g^^ is the nuclear splitting factor and = ^iTc" 
the nuclear Bohr magneton. Using first-ordor perturbation 
theory the eigen val'ues of "^(ml) are given by 

Sjj^(ml) = -ji H m^yj 

= ^ ( 2.10 ) 

This dipole interaction splits a nuclear state jl> into 21+1 
equally spaced states having nuclear magnetic spin quantum 
number m^ = I, The gamma transitions correspon- 

ding to these sublevels are subject to the selection rule 
1 = 1, Am = 0,+l. For ^"^Fe with I = ^ and = ■^ 

magnetic dipole interaction gives altogether six allowed 
transitions (Fig. 1.3). 

The relative intensities of the , six transitions are 


given hy 
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Fig 1-3 Magnetic hyperfinc splitting without and with electric qudrupole 
interaction . 
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( 2 . 11 ) 


For polycrystalline s.amples, the values shoiild be averaged over- 
all possible orientations of crystallites. This gives the 
intensity ratios to be 5c2;l;l;2:3. 

Pure nuclear raagretio dipole intora^ctions are rarely 
encountered in chemical a,pplications of the Mbsebauer effect. 

On the other hand the nuclear state is sinultaJ.eously acted upon 
by both magnetic dipole and electric quadrupole interactions. 

So the sublevels of I = ^ of ^ ^?e are no lo.uger equally spaced 
(Pig. 1 . 3 ). The qusdrupole coupling may be treated as a first- 
order perturbation on the- magnetic dipole interaction. If 
■^(el) ana if the mPO tensor is axially symmetric and 

its principal axis marcs “'n .-angle p with the axis of the 
magnetic field, from first-order perturbation calcul.ations one 
gets for the eigen values 

mj + ^ ? 

B = -gjj pj. Hm^ + k-l) (eQ (5Gos^p-l) 

( 2 . 12 ) 


1 . 3.5 Relaxation effects in Mbssbauor spoctroscopy 


When the environment of the Cossbauer nucleus rapidly 
changes, the M'dssbauer spectrum is affected. The fluctuations 
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may "be due to for example (ij time- dependent magnetic hyperfine 
interaction due to electronic spin-spin or spin-lattice relaxa- 
tion, ( ii ) randomly changing electric qup.drr.pole interaction due 
to Jahn-Teller distortions, (iii) interesti cial or vacancy 
diffusion in the vicinity of tho Mossbo,uer nucleus etc. 

Ihe lineshapo m.ay be quite- unueualbecause of rapidly 
changing Hamiltonian. Wo shall discuss in somewhat detail the 

effect of spin-lattice and vspin-spin relaxation on a quadrupole- 

^7 

magnetic spectrum of Fe. The proba.bili'Gy of emission of a 
photon in direction k with frequency e and polarisation p by a 
system which, in the process of emission, makes a. transition from 
its initial state I X | to the final state la> is given by the 


Wigner-Weisskopf formula [9] 


Xa 


(k 0) p) = 


<ai 7^ + ^ U>I 


( « + E,,. - Eyf -t- T^/4 


( 2 . 13 ) 


being the intcc cc^icr Eamiltonion be‘’‘V', in the nucleus end 
the el Gctromagnetiii field, j X > and la> represent in general, 
eigenstates of tho nucleus together with the entire solid or 
liqi^€’ in which the gamma-omitting nucleus is embedded. The 
quantity r is the inverse of the lifetime of the excited state 
I X> . Since we do not observe initial and final states in the 
usual experiments, the observed probability is obtained by 
summing over final states and averaging over iniitial states. 
Algebraic manipulations finally yield the anission probability 


[ 10 ] 
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* "P / \ 

P(ku)p)=-lR9 / d: [exp(iw t t)<7 

0 

.^fhere Re stands for real pa'’t, is the acRioii’t 

(t) = exp(i t) ' exp(-:. t) 


(t)>] 

of Y^'^^ and 


being the Haiiiltonian for ohe entire systtm. 

In a simple situation v:hero the electronic state of 

the ion containing the M'dssDauer nucleus flips between 

and - the nucleus finds itself in a magnetic field which is +H 

for sometime and -H for sometime. In addition, let us introduce 

« 

a constant axial electric field gradient, the principal axis 
being in the direction of +H, 

The time-dependent Haniiltonian may be written as 


-^(t) = ^ + f(t) + Qi3ll - I^) 


being the u 25 )erturbed Hamiltonian, 


o l"l“l" 


El 



E, 


1 m 
0 o 


> 


0 

and ^^“^o energy of the unperturbed gamma ray. The 

function f(t) can taXe two values +1 at random. The lineshape 
is obtained by averaging over all functions f(t). 


P(kajp) = •- Re / dr [expCloj (t)>] (2.14) 

’ o 

where the bar denotes stochastic average. The correlation 
function is calculated by inserting a complete set of states in 
the definition. 


\ 
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^T 


r( - ) tt( + ) / + 


(t) = 


21-,+l - . 

J. ^ Id. 

o . 


I <1.111. 


! V^"^l I m ><I m 
1 1 0 0 0 0 


21^1 V ^(t)ll^m^> 


where we have used the fact thst the photon emission operators can 
only connect the .excited, state I-,, ';‘ith tin. ground state 
Wri ting 


(t) 


exp(i f (t')c't’) exp(~.i / (t' )dt' ) 


t 


and inserting connlete sets i I ni' > <±ini i and jl.n’ > ! and 

° ^ '00 o' 'll 11' 

using 


<Iomo ! fexp(i / (t’)dt' ll^m;>= 6^^, oxp[i(B^+Q^{ 5mJ-I^(l^+l )>)t] 


0 0 


m m' 
0 0 


O O 0^0 


t 


X exp (ig^pEm^ / f(t’)dt’) 


and 


I^in£ exp(-i / ( t’ )dt' m‘ exr.[-i(B^+Q^{3m|-I^(l^+l)>) t] 

o ^ ^ 

t 


X exp (-ig^pHm^ / f(t')dt’) 


we get 


V^-) (t) exp (i AEt).-^j^ I m^>I^ 

1 m m. 

0 1 


X exp ( ip ( t ) )' (2.15) 

where Al = [3in^ - 1^(I^+1)] - E^~Q^[3in|-l3_(l^+l)] 

t 

p('t) = ^Sq^o - / f(t’) dt' 

o 

g^ and g^ are respectively the g-factors for the ground and 
excited states and and are proportional to the quadrupole 
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moments of the ground, and excited nuclear states. 

The matrix elements of v'' are proportional to Olehsch- 
Gordan coefficients and the only remaining problem is the 
evaluation of the average of the quantity 

JU 

0 

exp [i(gQm^ - g^m^)pH / f(t’ ) dt’ ] 


The average of this quantity may he expressed in tei*ms of the 
proha'bilities V/ arid V/ , of transitions from +1/2 to -1/2 and 
from -1/2 to +1/2 respectively. For = V/ 


exp[i(g m -g m., )pH / f(t')dt' = [(Cos x}Ut +“(Sin x)Wt] exp(-Wt ) 




X 


( 2 . 16 ) 


/ . 2rr2 

where x = [ 5 IJ 


Substituting eqs. (2,15) end (2.16) in (2.14) '"ives the line- 


shape. 

Spectra v:ere calculo.ted [11 ] using the above theory for 
different values of she jump time t. For large x or slow 
jiimping the spectrum corresponds to a fiill six line one as if it 
is in steady magnetic field. For very small t, the spectrim is 
a symmetric quadrupole doublet. For intermediate jirap times we 
get many lineshapes. 

As we start from large relaxation time x, full six line spec- 
trum corresponding to static magnetic field H with quadrupole 
modification is observed. As x becomes shorter the inner magnetic 
lines collapse onto their centre of gravity. The other magnetic 
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lines axe still present, altiicngh they axe considerably broadened. 

The renaining magnetic lines also go on collapsing and quadrupole 

Pc^ttern starts to snpear as x becomes shorter and shorter, [phe 
3 1 

"rt-line i.e. + ^ more slowly than the cr~line 

1 1 . 1 

±2 i ‘2* ± ■f ±2 "transitions, v;hich make up the m-line 

of quadrupole dou'Dlet, have a larger splitting in a magnetic 
field than the + ^ ^ transitions. The Zeeman precession 

frequencies, with which relaxation frequency is compai'ed rj:'e 
different for different transitions. It is possible that the 
relaxation time t is such that the + ^ ^ transition feels 

a rapid relaxation but + -^ ± f feels a comparatively slow 

relaxation. The high velocity components will have then 
already collapsed to a sharpi line whereas the low velocity compo- 
nents will still be broad. This then results in asymmetric 
quadrupole doublet, hecem^ly, c,aiculationG v-ere reporiud [12-13] 
where the effecic of "he motion of the Mo'^^s'*' 'U er nucleus itself 
or that of the mole-u] es (in case of liquids) etc. have been 
takien into account. 

1.3.6 Goldnski-Karyagin effect 

The intensity of the M'dssbauer line depends on the 
recoilfree fraction 

f = exp[- <(k.x)^>] 

* 

If the gamma ray direction k is expressed in polar co-ordinates Q 
and 0 with respect to EFG" axes, we can write 
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{ [cx^ > Oos ‘'’^0 + <z^> Sin‘^ 0 ] Sia^© + <x^> Gos^G ] 

If tile amplitudes of vibration are not equal in all the 
directions, the recoilfree fraction becomes a function of 
( 9 , 0 ), For ^ ^Fe quadrupole doublet for a polycrj^-s talline 
sample the intensities of the t\fo lines have tiie ratio 

2ti X P - 

/ / f(©, 0 ) ( 1 +Cos^©) Sin© d© d 0 

^ _ o 'o 

]}/[ 2 !.% 

^ I / f(©, 0 ) (| - Gos^©)Sin© d© d 0 

0 0 ^ 

for axially symmetric vibrations where 


^2 

‘z 


<Z„ > = <X^ > = < X J > :j=. < > =: < x‘i I > 


^ i| 

f( 0 , 0 ) = exp [-k^<x? >] exp[k^( <z‘T> - <;-:?■ > ],Oos^©] 


resulting in 


/ 3npr-l;^(<mp >- <x^ > )Gos‘'©] vl+Oos^© )Sin© d© 

% o [ 

^ I exp[-k^( <xt. > - <z^ > )0os‘"'|(’4 - Cos^©)Sin© d© 


We thus see thp.t evpoi for random orientation of crysta- 
llites ,the intensities of the two lines ma^y'" be different. The 
intensities of the magnetically split linos arc also modified - 
for the same reasons, i.e. the anisotropy of the De by e-Waller 
factor. This effect is known as Goldanski-Karyagin (G-K) effect. 
As the vibrational anisoti'opy is an increasing llinction of 
temperature the asymmetry in intensities related to the G-K 
effect also increases with temperature. 
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1.3.7 Asyiometric quadrupole doublet 

The observed quadrupole doublet may be asyinmetric due 
to the following three reasorss 

1) lack of randomness of crystallites 

2) Anisotropy of Debye-Wallei’ factor 

3) Relaxation effects. 

The first one has the distinct propt.rfcy that the ratio 
of intensities strongly depends upon the angle between the 
gamma ray direction and the sample surface. The second and 
third effects are independent of sample drientation but depends 
very much on temperatxire. The intensity ratio due to anisotropy 
of the Debye-Wall or factor is ^1 at low temperatures and 
increases with temperature. For the third case the asymmetry 
ratio will strongly depend upon the relaxation mechanism. If 
the fluctuations in the magnetic field are due to spin-lattice 
relaxation we expe-r-t relaxation time to decrease with increasing 
temperature. This me^^ns tli..t raising the temperature will 
decrease the asymmetry, an cilect opposite to temperature depen- 
dence of G-K eifect. On tie other hand if the fluctuations 
are due to spin-spin rela:cation, the asymmetry will depend on 
the param^netic ions. 
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CHAPTER 2 


EXPERIMENTAL 


The principal components and equipment used in these 
experiments were ; (i) a source of y-raj generated by the decay 
of 'Go to ^'Pe,(ii) the sample containing ^ ¥e, always used 
as the absorber, (iii) an electromechanical system acting as a 
drive unit to impart velocity to the source, (iv) the dewar 
and the furnace for keeping the absorber at the required tempe- 
rature, (v) a detector and (vi) a storage device (the 

multichannel analyzer) which stores the data as a function of 
the velocity of the source. 


2,1 Instruments 
2 ,1 .1 Electronic equipment 

a) Mbssbauer spectrometer 

The spectrometer used in the present work contains a 
programmable velocity drive and an electromechanical transducer 
which was developed in this laboratory and constructed by M/s 
Encardio-Rit e-Lucknow, This spectrometer has provision for 
constant velocity, constant acceleration and sinusoidal velocity 
modes. A precision velocity drive is realised, by generating 
signals with highly stable amplitude, frequency and phase, and 
connecting the transducer in a feed back loop. The block 
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diagram of the circuit is shown in Jig, 2.1. The drive signal 
is first genei-ated in a higiily stable 100 kHz crystal controlled 
oscillator. The output of this oscillator is scaled down by 
JK Plip-ELop to a 20 Hz square wave. The reference triangular 
wave for constant acceleration mode is generated by integrating 
the square wave by an active integrator using 10 741, The 
reference sine wave for the sinusoidal mode of the drive is 
obtained by extracting the fundamental component of the square 
wave through an active filter. The use of square wave to the 
input of the active filter minimises the second harmonic distor- 
tion in the output .because of the absence of even harmonics in 
a square wave. A reference trapezoidal wave for the constant 
velocity mode is generated by integrating the square wave in a 
suitable integrator. A trapezoidal wave, instead of square 
Wave is used, so that the slew rate of the power amplifier and 
transducer is not exceeded which would otherwise result in 
ringing. An automatic gain control (AGG) circuit is provided 
after the integrator/ active filter to stabilize the amplitude 
which is very critical for good performance. The AGG output is 
used as the reference signal for the transducer drive coil. 

The transducer consists of a conventional double speaker 
- arrangement .having a drive coil and a pick-up coil. The velocity 
is monitored by comparing the pick-up signal with the reference 
signal in an error amplifier. The difference between the two 
signals, caUed the error signal, is a measure of the deviation 
from the desired velocity. The error amplifier output is fed 
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to a power amplifier which, drives the transducer, Thus the 
transducer velocity follows the reference signal precisely in 
amplitude and phase. An MCS start pulse to reset the sweep in 
the MCA to the first channel is obtained by scaling down the 
complementary output of the penultimate -7 2 stage by a factor 
of 2 which is again applied to a monostable multivibrator 
(74121) through a level shifter to get a narrow pulse of width 
10 |j.sec. The use of the complementary output synchronises the 
opening of the 1st channel of' the M0A with the negative maxima 
of the vej.pcity reference signal. 

The design of the transducer is shown in Pig, 2,2. The 
impedance of the drive coil is matched with the output impedance 
of the power amplifier. The transducer is so designed that the 
frequency of unwanted vibrations, e.g, drum vibrations or 
overtones of the guided springs etc, are far away from the basic 
resonance frequency. The pick-up coil and drive coil are wound 
in grooves, at the two ends of a common holder, having diameter 
of 1.5' which is made out of perspex. Pach coil has a length 
of 0.3’’ which is subdivided into three coils of length 0.1’’ 
each. The three pick-up coils are connected in series and the 
drive coils are in parallel in the normal operation for impedance 
matching of the drive coil and to improve the linearity. 

B* Detector and amiplifier 

Por detection of y-rays in our experimmts we have used 
an argon-filled proportional counter with alun^nised Mylar 
Window, procured from the Technical Physics Division, BURG, 



Source Mount Permanent Ring magnet 



chcmatic diagram of the electromechanical transducer 
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Trombay, Bombay and a .Krypton- filled proportional coimter 
man^ifactured by Electronic Corporation of India ltd. The 
author has fabricated a low noise, highly stable , temperature 
compensated ,pre- amplifier for the proportional counter the 
schematic diagram of which is shown in Pig. 2.3. The signal 
from the pre-amplifier was amplified by an ORTBC Model 410 
linear amplifier, the output of which is fed to the ED-60 
multichannel analyzer. The P¥HM observed for the 14 keV gamma 
rays with this set up is 2 keV for the argon- filled propor- 
tional counter. 

C. Multichannel analyzer (MOa) 

The NI)-60 microprocessor based 2048 channel miilti- 
channel analyzer was used for selecting the 14.4 keV gamma 
rays and the storage of data. The MCA is operated in the 
multi-scaling mode with external channel reset, with a dwell 
time per channel of 100 psec for 512 channels operation and 
50 psec for 1024 channel storage. 

Multichannel scaling is comparable to counting in a 
series of scalers each being put into operation successively 
for a preset time, called the dwell time per channel. The 
address register in ED— 60 can be advanced by the internal 
time base with a minimum loss of counts por channel during 
transition, so that the effect ency of coimting is much higher, 

f 

The sweep is triggered' by the start pulse f 
in the MGS external mode, to get a stable p; 
acquired in each subsequent sweep is added ■ 
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stored in tiie previous sweep. 

The use of microprocessor memory in this analyser 
gives a better counting efficiency compared to the magnetic 
memory storage. The data can be printed out in a series line 
printer (Teletype) or can be plotted in a Hewlett Packard 
Model 7015A X-Y recorder, 

D. Performance 

The performance of the system was found to be good. 

The full width at half-maximum of the Mbssbauer line obtained 

in the present work: with Pottasium Perrocynide absorber and 
57 

Co in Rh matrix source was 0.27 mm/sec. This width 
includes the contributions from thickness of the source and 
absorber etc. The velocity per channel was found to be 
reproducible with in 0.1 over a period of several weeks, 
once the spectrometer is set for a particular calibration, 

2.1.2 Furnace 

The designs of vacuum furnaces for Mbssbauer work at 
high temperatures ( '^ 300°-1300°E) are available in the lite- 
rature [1,2, 5,4], Usually the heating coil is made of kanthal 
D wire wound on an alumina muffle mounted by means of ceramic 
flanges inside a brass cylinder (about 15 cm long and 15 cm 
diameter). Because of the involvement of large mass and 
different materials inside the furnace a large amount of power 
('w 1 KW) is required for these furnaces and the amount of 
degassing and contamination of the sample is large at higher 
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temperatures. The outer brass cylinder requires an elaborate 
Water cooling arrangement. The large size of the furnace brings 
restrictions on the versatility of its use, e.g, the geometry 
becomes poor with week sources, and giving to and fro motions 
to the sample for recoiless fraction measurements etc. becomes 
difficult . 

The author has designed a miniature quartz vacuum 
furnace [5] of light weight, requiring only small power 
{ 40 Watts), with temperature controlling facilities, and a 
small size giving good geometry and offering flexibility for 
various experimental operations. A special feature of the 
furnace is that it has got arrangements for induction as well as 
resistance heating. The induction method reduces the contami- 
nation of the sample, minimises heat losses and gives better 
temperature stability throughout the range of operation. The 
two modes of operation, i.e. induction and resistance, can be 
very easily interchanged , The resistance method may be used 
upto lOOO’^K whereas the induction method raise the sample 
temperature to 1300°K and, with suitable modifications, even 
higher , 

The design of the furnace for resistance heating is 
shown in iig, 2,4. A quartz tube of diameter 4,5 cm and 
length 8 cm has one end tapered to a diameter 2 cm. This end 
is permanently sealed with a mylar window which is shielded 

from thermal radiations of the heating coil by a thin sheet 

\ 

of beryllium, A flange is made at the other end, on which a 
circular brass or perspex disc with a mylar window (of diameter 





2 cm) at its center, is kept pressed with an 0-ring in 
between. The thermocouple and the heater connections are 
taken through the brass/perspex disc, A copper tube of diameter 
6 mm is wound around the quartz tube for use as water-cooling 
system. This coil is also used as work coil during the 
induction heating. The tube can be evacuated by connecting 
a vacuum pump through the side tube, which can be closed by 
the stop cock after evacuation. 

The heating assembly consists of a quartz or alumina 

muffle of 3 cm diameter and length 4 cm with nichr'ome wire 

wound on it non-inductively . The coil is kept in position 

by high temperature cement. Spacers made out of refractory 

brick are fixed at the two ends of the muffle. 

\ 

The sample holder in the case of resistance heating 
is a copper disc of diameter 2.3 cm and thickness 1 cm with 
a groove of 2 cm diameter and an axial hole of 1 cm diameter. 

The M’ossbauer absorber is kept between two thin discs of Pe- 
free synthetic mica, put on the groove and pressed by a copper 
ring. The sample holder is pushed snugly at the center of the 
muffle. The entire heating Glomv.-nt with the sample is now 
put inside the outer quartz tube and the perspex/brass side 
disc clamped with brass clips. The temperature of the absorber 
is sensed by a chromel-alumel thermocouple, and kept steady 
by a solid state temperature-controller. 

For induction heating the resistance heating system 
inside the quartz tube is replaced by an alumina tube of 
diameter 3 cm kept in position by refractoiy brick spacers. 
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The saiaple holder is made of two tantalum discs of diameter 
2,5 ^m and thickness 1 mm, kept in position inside the alumina 
tube by alumina rings. The high frequency (200 kHz) ciorrent 
is passed through the outer copper coil generating an induced 
eddy current in the tantalum sample holder and thus Joule 
heating the sample. The temperature is controlled by the 
temperature^-controller which can switch on or off the high 
frequency current. The induction heating power supply is 
basically a plate- tuned, vacuum tube (807 tetrodes) feed-back 
oscillator with the work coil coming as a part of the inductor 
in the tank circuit [Fig. 2.5]. 

The temperature-controller which the author has fabri- 
cated is basically a chopper stabilized amplifier working on a 
proportional on-off principle with the circuit as shown in 
Figs. 2,6a and 2.6b. When stabilized the load is switched 
on and off once every second. The duty cycle is controlled by 
the difference between the thermocuple and the reference levels. 
As the response time of the furnace is very long the furnace 
itself acts as an integrating element and the observed tempera- 
ture is due to the average power being put to the system. For 
a miniature furnace as described above the temperature can be 
regulated to 0.2°G for low temperatures and within 0,5^0 for 
high temperature. The characteristic of the design of the 
controller is such that it can be used with a load of any 
impedance. 

The furnace systaa with the accessories has been used 
extensively in our laboratory for the studies of the toaperatur© 
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vai'iation of Mbssbauer spectra of a variety of systems. For 
the heating coil arrangement a 5 amp. -50 V d.c. commercial 
power supply was used. For a 50 Watt supply a temperature of 
700°G could be attained. For the induction furnace an RP power 
of 250 Watts raised the temperature to 900^0. The temperature 
can be further raised by keeping the work coil inside the 
furnace. The temperature-controller was calibrated with 
respect to a chromel-alumel thermocouple by adjusting the refere- 
nce voltage in the circuit so that each division of the ten- 
turn potentiometer dial corresponds to 1°G . 

Some typical examples of the operation of the furnace 
in our laboratory are ; study of the magnetic anisotropy and 
thermal decomposition of Fe(Rh:^)2(S0^)2.6H20 single crystals 
( 300®K-500®K, with heating coil), study of Mq 57 alloys 

at different temperatures ( 500°K-500°K, with heating coil), 
study of Fe-Pd alloys ( 300°K-900°E, with induction heating) etc. 
It is to be noted that there is no sideband by RF excitation 
[6], This is due to the fact that the frequency used is very 
small (*'' 200 kHz) compared to that required 30 MHz) for RF 
excitation of magnetic materials [7]. 

2.1.3 Gryostat 

A liquid nitrogen dewar was fabricated by the author for 
taking the Mhasbauer spectra at t^perature between 80-300 K. 

The design is shown in Pig. 2.7. It consists of two concentric 
cylindrical brass vessels both of which are connected by a thin 
walled stainless steel tube through which liquid nitrogen is 
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poured into the inner container. A demountable copper cold 
finger is attached to the inner vessel with provision for 
introducing suitable thermal insvilation material between them, 
for efficient stabilisation of temperature of the cold finger 
at temperatures higher than 80°K. The thermal insulation 
between the two containers is provided by evacuating the space 
between them to a vacuum .02 micron which is measured by the 
ionisation gauge connected to the dewar at the top. The 
absorber can be kept pressed to the cold finger by flat thick 
copper square plate. This gives a better uniformity of tempera- 
ture over the sample zone. Plat nichromo wire is .wotiodon the 
cold finger with mica insiuLation near the sample zone for 
heating the finger. The temperature of the sample is measured 
by means of a copper-constantan thermocouple. The leads of 
the thermocouple and that of the heater are taken out through 
glass to metal seals fixed to the outer vessel. Demountable 
windows are provided on the outer vessel assembly for admi- 
tting gamma rays. Because of the demount ability any type of 
window for the particular experiment can be easily fixed. 

The bottom flange can be easily opened to change the 
samples. All parts of the dewar are nickel plated to prevent 
heating by radiation. A stainless steel stand assembly with 
rubber mountings for the legs^-S designed, and the dewar can be 
fixed to it at any height by the screws provided. 

The temperature controller described along with the 
furnace is used to keep the temperature constant at any value 
between 80-500 within an accuracy of 0,4 °C. With one 
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poured into the inner container, A demountable copper cold 
finger is attached to the inner vessel with provision for 
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the thenno couple and that of the heater are taken out through 
glass to metal seals fixed to the outer vessel. Demountable 
windows are provided on the outer vessel assembly for admi- 
tting gamma rays. Because of the demount ability any type of 
window for the particular experiment can be easily fixed. 

The bottom flange can be easily opened to change the 
samples. All parts of the dewar are nickel plated to prevent 
heating by radiation. A stainless steel stand assembly with 
rubber mountings for the leg^s designed, and the dewar can be 
fixed to it at any height by the screws provided. 

The temperature controller described along with the 
furnace is used to keep the temperature constant at any value 
between 80^300 within an accuracy of 0,4 ^0, With one 
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filling of the dewar by liquid nitrogen the liquid remains for 
about 12 to 14 hours. 


2.2 Sources and Absorbers 


2,2®1 Sources 
57 

The Co source in Rh matrix of strength 20 mCi used 

in these experiments was obtained commercially from Uew England 

Nuclear, Massachusetts USA. This source gave a line width of 

0,26 mm/sec for Pottaaium Perrocynide absorber >ath 0,5 mg 
57 ! 2 

Fe/cm distribution, with our set up. This width includes 
the width due to source, absorber and broadening due to absorber 
thickness. The recoilless fraction of the 14.4 keV y-^^ays from 
this source is found to be 0,78 »as measured by Trooster et al. 
[8] method, using natural iron foil 0.0005” thick and a brass 
filter 0.005” thick [9] for background correction. 

This source satisfied the following requirements: 

i) The host matrix is chemically stable. 

ii ) There is no electric field gradient or magnetic field 
at the site of the Mbssbauor nuclei, 

iii ) The recoil- free fraction is large, 

iv) The line width introduced by self absorption, thickness 
effect etc. is fairly small. 

v) The charactOTistic X-rays from the source matrix Rhodim 
(20-23.3 ke7) is well above the Mbssbauer gamma radiation, 
so that they can be easily eliminated by the single channel 
analyser. 



t 


(2,3) 


= CfH f 

a 

where a is resonance cross section of the absorbing nuclei 

S 7 V 2 

(^^Fe), n is number of absorbing nuclei/ cm in the path of 

Y-rays and f^ is the recoilless fraction of the absorber. 

In order that the least time is spent in M'ossbauer 

absorption measurements, an absorber should be chosen, so that 

it will provide optimum statistics in a given time. That is it 

increases resonance absorption and reduces non-resonance 

absorption as far as possible. 

The number of transmitted gamma rays at the resonance 

dip as a function of absorber thickness is given by the 

equation [9] 

N(t,V^) = (1-f) N + m exp(-t/2) J^( i t/2 ) (2.4) 

where t is the thickness parameter, is the Doppler velocity 
and N is the number of Y~2*ays incident corresponding to the 
particular transition in the absorber. The amplitude of the 
Mbssbauer line is then 

= m [l-erp(-t/2) (2.5) 

J^(x) being the zero order Bessel function of imaginary 
argument. The above equation expresses only the resonance 
effect as a function of Doppler energy. This expression 
suggests that the resonance absorption amplitude increases 
with thickness mo no tonic ally . It is obvious that as the 
number of active nuclei in the path of the Y-^ay photon 
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increases the absorption will be more. However, if the 
absorber, is too thick the non-resonant damping of all gajama 
rays will increase and the background noise will also increase 
resulting in masking of the spectrum. On the other hand 
making the absorber too thin, although increases the emergent 
beam, would reduce the resonance effect because of the fewer 
active nuclei in the path. Thus there must be an intermediate 
thickness, which for a given number of incident photons, will 
give a large absorption peak. When the non-rosonant absorp- 
tion is also taken into account [13] the amplitude becomes 

N-H{d,V^) = fN[l-Gxp(-^ij,d/2) J^(i ■^)] expC-Pg^d) 

( 2 . 6 ) 

when u and are the mazimina resonant and non-resonant 
absorption coefficients. The relation between and o is 

l^r = (f) ^ (2.7) 

when d is the actual thickness of the sample and ^ is the 
density of the active nuclei. Shimony [14] has maximised 
this with respect to d to get the highest peak for a given H. 
The optimum thickness in that case is given by the relation . 

Iij. ezpd/a d^)-j^(i/2 i 

Jq( 1/2 i d^) +Jj_{l/2 i d^) 

d(J.(z)) 

using the identify = -J^(z), 

If we instead concentrate on a good spectrum with 
minipum statistical uncertainty a better parameter for that is 
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Amplitude of the absorption dip 


Statistical uncertainty in the counts 


At velocities far away from the counts are 
N(d, “ ) ?= N e2?;p(-.u d) 

S- 


If the absorption is not too larg-e the statistical uncertainty 
is ^iven as 

= [N(d, 

= [N exp(T.pg^d)]'2 

, , a(d) = f \/'N' [3,^exp(-[i^d/2)J^(-r-^^)] exp(-p^d/2) (2.9) 


Equating the derivative of a(d) with respect to d to zero we 
get 


p exp (1/2 pj. d^)~J^(l/2 i p^ d ) 

n^-'^WTOWi.^Triitcnn-t- - ii — Ti i Ti HI i ii «iii mp rhiMiirrri^a.irf|-riiiTlrWtMirihii.Htir irii[-nh-n..T.L- ni.mi-r it— . -rr .-ii^wirr 

Ha ■■■ Jq( 1/2 i p^ d^)+J3_(l/2 i Pr 


( 2 . 10 ) 


from this equation d can be found out in terms of p^, and p . 

U X 3- 

In the case of EeS0^,7H20 using the above method and 
Pj, = 35.96 and p^^ = 22 <ind from • the table in ref. [14], -we ^t 

d ^ 40 mg/ cm^. 


2,3 Effect of G-eometry 

< 4 

The geometry introduces two types of distortions in the 
M^di^sbauer spectrum. Because of the movement of the source the 
distance between the source and detector changes. The cotint 
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rate then changes with an 1/^ dependence which causes a non- 
flat base line in the spectriun. For exajaple in the constant 
acceleration mode the velocity at any instant t is given by 

V(t) = ^ t 


for 0 t <.T/2 and 

V(t) = " -T|^t for T/2 <. t 

where is the amplitude of velocity and its time period. 

The distance between source and absorber at any instant from these 
equations is 

0 < t <. T/ 2 and 
T/2 <.t j< T 

Xq is the distance of the source from detector at t = 0, i,e, 
when the source is at its equilibrium position^ At t = T/2 


x(t) = -V^t + t^ + 

27 

x(t) = i^^pt T t^ + Xq 


X 


o 


o o 


V/riting x(t) = + x(t) 


x(t) = 

-v^t + 

2^0 0 
0 f2 

■T' ^ 

vl 

0 

t 1 t/2 

and 

xCt) = 

2V„t - 

2V ^ 

0 4.2 

"T“ ^ 

T/2 

1 t < T 



The count rate in the absence of any absorber is given by 


G(t) 


K 

[Xp + '^X( t 


where K is a constant 
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The maximum value of |az{ is at t = T/4 or 3T/4s"‘i'hero the 
velocity is zero, i,e, at xhe turning poinuSo 


K 


ax)! 


max 


V T 
0 

S 


In our system for T = 50 millisec. and ^ - 

0.05 3am. If x = 5 cm, ™ is nearly 0<,0''l. 

^o 

expression for C(t) using hinomial theore ' 


V T 

10 mm/sec, -g- = 
By expanding the 


0 ( t ) = [ 1 — 

X ^ ^0 

o 

The first order shift in base line then comes as 0.1 / . , 

This Can be neglected for samples having good percentage 
absorption.. In our experiments we have kept the distance 
between the source and detector at 10 cm. The method of 
nullifying this effect is to ’fold* the spectrua about the 
maximum velocity. This eliminates the effects due to ohange 
in count rate in the fir.et coder and impro’-es, the statistics 
:.f the lines are nsirrow l- 15]- lines.rity of 

Ghe Velocity cycle is i: orooed in bo -u x.aa.v6S by this 
method.^ " ^ . 

The other effect known as the ’cosine smearing effect* 
is caused by the finite size of the source and the absorber. 

The gamma rays reaching the detector have different angles 
with respect to the direction of motion of the source. If this 
angle is ©, the gamma ray energy as seen by the absorber will 
be 

^ E (1 . nSaafi) 

O Q ^ 




E’ (e) 
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If 6 is large (when the detector and source separation is small) 
this effect will be pronounced and it causes a broadening of 
the lines in the absorption spectr’am and even changes the shape 
of the spectrum. 

2„4 Data Analysis 

2.4.1 Oalibration 

Por finding the velocity per channel of the sample, 
calibration spectra were tahen before and after each experi- 
mental run, lor calibration samples were either enriched iron 
or sodiiim nitroprusside (SNP) in our experiments. 

For isomer shift measurements SIIP and 310 stainless 
steel were used as reference materials. These standard 
absorbers were obtained commorcially from Radio Chemical Gohtro 
Amersham, Harwell, U.K. 

2.4.2 Computer programmes 

If thu relaxation effects are small, the Mbssabuca* 
absorption line can be well approximated by a Lorentzian curve. 
There are many programmes available now to analyse the 
Mbssbauer data. The basic method is almost the same, where 
some functional form for the absorption peak is assumed, frr 
example as lorentzian or Lorentzian + Gaussian etc. The man 
purpose of these curves fittings are to find the values of t'^'e 
adjustable parameters (g»Aess values), which provide the clcsv^it 
agreement between experimental data and the values predicted 
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by the assumed function. 


The method us :d to find the best panpireter in our 
programme is based upon the method developed by V.J. law and 
R.V, Bailey [16,17]. Let 0^ ■= be the 

functional form chosen for the experimental points Y. We 
wish to find a set of parameters such that 


N 

I 

i=l 


S. 




I Ui 


i=l 


0 ,) 


c. 


is minised (N ie the number of data points). The condition 
for this is (dS /Bj^) r= 0, If this is done a set of m simul- 
taneous equations will be obtained „ These equations are 
solved tp get the corrections aBj^. and new values of the 
parameters obtained where tho super- 

V"' 

script i refers to the i iteration. 

The ensure that the new values are better, AB-^ is 

2 

chamged in small steps at the point of convergence so that S 


is minimuBi. 

The function 0 and its derivatives ai’e calculated in 
a separate subroutine, 1 chatige in function, al form can be 
accomodated by making changos only in this subroutine. 

The progrojame also gives the error matrix [18 j from 
which the standard error in the fitted parameters B„ can be 
calculatede All the computation are done in the DEG 1090 
computer of I.I.T. Kanpur, We have developed other computer 
programmes for finding , the P(H) distribution in magnetic 
hyper fine split spectra as well as BPS distribution 
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in the case of quadrupole split spectra, which will he 

discussed later in this thesis* Secoilloss fraction and Dehye 
temperature of some of the samples were also found out for 
which computer programmes were made hy the author. 



REFERENCES 


53 


[1] D.E, Nagle, Phys. Rev. Lett. J, 564 (i960). 

[2] R.S. Preston, S.S. Hanna and J. Heberle, Phys. Rev. 128 . 
2207 (1964). 

[3] V.S, Sundaram and V.P. Gupta, Rev. Sci, Instr. _42, I6l6 
(1971). 

[4] F. Van der Woude and G. Boom, Rev. Sci, Instr. 36, 800 
(1965). 

[5] N.V. Nair and D.G, Khan, Nucl. Instr. and Meth. ( comm- 
nicated) . 

[6] B, Pfeiffer and N.D. Heiman, Phys. Rev. m, 74 (1972). 

[7] M. Chmielowshi and A. Kotlicki, Hyp, Int. 1129 (1981). 

[8] J.M, Trooster and M.P.A, Viegers, ’’Mossbauer Effect 
Reference and Data Journal” ,1, 154 (1978), edited by 
J.G. Steven-, (I FI /Plenum, New York). 

[9] R.M, Housley, N.E, Erickson and J.G, Dasn, Nucl, Instr. 
and Meth. 29 (1964). 

[10] S, Margulies and J.R. Ehrman, Nucl, Instr. and Meth, 12 , 
131 (1961). 

[11] D.A. Shirley, M. Kaplan and P, Axel, Phys, Rev. 123, 816 
(1961). 

[12] D.A. O'Connor, Nucl. Instr. and Meth. 318 (1963)* 

[13] A,J. Stone, Nucl. Instr. and Meth, 107 . 285 (1973). 

[14] U. Shimony , Nucl. Instr, and Meth. 348 (1965). 

[15] lakayiiki Kobayashi' ,and Kaztiko Fufcumura, Nucl, Instr. 
and Meth. l^^, 363 (1980). 



54 


[16] VeJ. Law and E..V, Bailey, Chem, Lng, Sc„ 1^, in9 (1963). 

[17] H.O. Verma, Ph.'„D<, Thesis, Indian Institute of Technology, 
Jfanpur, India (1980) (unpublished,, 

[IS] J.R. Wolber-g, ‘ predioti on Analysis ’ , L, Van Bostrand 
Oompany Inc^, Princeton (1967). 



GHIPTER 3 


MOSSBAUER STUDIES OE SOME FERROUS COMPLEXES - A RE¥ 

LOOK lETO THE OLD SYSTEMS 

3.1 Introduction 

We present in this chapter some new findings in the 
well known systems such as PeS0^.7H20, Ee( 30 ^ )2. )2* 6H2O , 

when an electric field is applied parallel to the direction 
of [1] and also the electric field gradient distribu- 

tion and isomer shift distribution in several such complexes 
using the model independent method. The absolute recoilless 
fraction in these systems are determined using the method 
adopted by Johnson et al, [ 2 ], with mo fifi cations. 

3.2 Mhssbauer Studies of PeS0^.7H20 and Fe(S04)2(MI^ )2. 6H2O 
Under Applied Electric Field 

Mhssbauer spectra of powdered samples of PeS0^.7H20 and 
Fe( 30 ^ )2(UH^ )2» 6^2© under an externally applied field along 
the direction of the were studied. Even though the 

internal crystal field was expected to be several orders of 
magnitude more than the applied field, the detection of the 
effect by a high-resolution spectrometer was not ruled out. 

As reported below, the indtjced St^ll effect was clearly 
identifiable in one of the systems, namely, FeS0^,7H20. 



56 


Commercially available high-purity FeS0^,7H20 was powdered 

(Merck analytical reagent) and kept between two capacitor plates 

2 

as a thin layer with a distribution of 40 mg/cm . Mbssbauer 

spectra were obtained at room temperature for various applied 

fields up to 45 kY/cm along the y-ray direction using the 

highly stable constant-acceleration spectrometer. 

The Mbssbauer absorption spectr-um of PeS0^.7H20 at zero 

external field was a quadrupole doublet with AEq = 3.198+0.003 

mm/sec and isomer shift 1.383 + 0,002 mm/sec [Pig. 3.1( a) ] . The 

values agree quite well with the standard values quoted in the 

literature [3-7]. Under the application of the external 

electric field, it was found that AEq increases up to 

3.247 + 0.003 for a maximum field of 45 kV/cm [Pig. 3.1(b)], 

The variation of AEq with the applied field is shown in Pig. 

3.2, The line shape is' Lorentzian with full width at half 

maximum (P¥HM) 0.278 + 0.004 mm/sec (for zero field) which is 

57 

comparable with linewidth of Co in Rh. 

The isomer shift changes from 1.383 + 0.002 mm/sec to 
1.411 + 0.002 mm/sec for the same range of applied voltage 
(Pig, 3.3). The intensity ratio of the two peaks (l-^j + 2 ^ 

1-^, + ■g> and + ■^> transitions) at room tempera- 

ture at zero field was unity. However, with increase of field 
a texture effect was produced which generated an asymmetry of 
line shapes. The intensity ratio decreased linearly with the 
field, reaching a value of 0.86 at 45 kV/cm (Pig, 3.4). The 
texture effect cannot be due to alignment of the small crystallites 
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FigSl Mossbouer spectra of FeS 047 H 20 at 
fa) Zero field (b) 45KV/cm. 
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Fig. 3-2. Change In quadrupole splitting with 

applied field In Fe 50^ THjO (curve a) 
and In Fe(504)i (NH^), SHjO (curve b) 
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Fig. 3-4 Variation of inttnsHy ratio of fht 
two pecks with opplied field In 
Ft S0^.7Hj0 



of the powder, since they were immobilized with an epoxy binder. 
The effects described in the Pigs. 3.2 3.3 and 3.4 were found to 

be reversible as the external electric field was decreased from 


V to zero, 
max 


The experimental values at room temperature for 
Fe(NH^ )2( SO^ )2 . 6H2O at zero field were 1.764 + 0.002 mm/ sec 


for quadrupole splitting and 1.400 + 0.002 mm/sec for isomer 

shift in agreement with the standard values [6,8]. However, 

% 

there was no change in any of the above values up to an applied 
field of 31.4 kY/cm, Above this field two new outer peaks 
appeared with AEq = 2.885 + 0.002 mm/sec, and isomer shift 
1.317 + 0,002 mm/sec. The spectrum for this system at zero and 
31.4 kV/cm are shown in Pigs, 3.5(a) and 3.5(b), respectively, 

A texture effect was observed for this system also, 

Bashkirov and Lebedev [9] have recently made a pertur- 
bation calculation of the polarization of the Pe^’*’ ion 4s shell 
by an electric f ield and its contribution to the quadrupole 
splitting and isomer shift, Por a field E = 45 x 10^ Y/cm, 
using their Bqs, (2) and (3)> the 4s-shell polarization effect 
tends to zero. Hence, the contribution could be due to the 3d- 
shell polarization and/or the effect of the induced dipoles at 
the lattice sites (for aEq only). To our knowledge, the 
problem of the effect of a perturbation Ez on the 3d-electron 
states has not been worked out yet. Another approach would 
be to use the axial symmetry of the operator and suitably 
augment the coefficient of Eq. (7) of Iiigalls [3] paper on 
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Fig3S. Mossbauer spectra of Fe(S 04 )i{NH 4 ) 2 ‘ 6 H 20 
at (a) Zero field (b) 30KV/cm. 
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EFG tensor in ferrous compounds. The lattice dipole effect 
could be calculated classically. Experiments with other 
systems with larger fields up to the dielectric breakdown 
limit will enhance the requirement for a detailed theoretical 
calculation. 

The extra peaks (Fig. 3.6) in the Fe(S0^)2(KH4)2‘^H2*^ 
appearing with the application of high voltage sufficient to 
break the dielectric is possibly linked either with the change 
of the hydration states [10] or v/ith the implantation of 
electrode (Al) atoms in the system. 

3.3 Electric Field Gradient and Isomer Shift 
Distributions from Mbssbauer Spectra 

3.3.1 Introduction 

Recently considerable attention has been paid to the 
line broadening of the M’dssbauer spectra due to absorber 
thickness [11-22], But even with very thin absorbers the line 
broadening can be appreciable due to the distribution of the 
isomer shift, the quadrupole splitting and the hyperfine field. 
Chien [23] has applied Window* s method [24] for finding the 
hyperfine magnetic field distribution in amorphous ferromagnet 
FegQB 2 Q and Rarayanaswamy et al. [ 25 ] applied this method to 
find out hyperfine field distribution in disordered binary 
alloys of Ni^Fe, RiFe and FeCo. Eibschutz et al. [26,27] 
interpreted the large width of vitreons YIG as being due to the 
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distribution of the electric field gradient at Fe sites. They 
described the quadrupole split doublet line shape f(Z) by a 
distribution p(lvl) of lorentzian lines. The EFG distribution 
p(lVl) Was derived from a computer model developed for the 
sequential addition of hard spheres [28], 

The present authors have applied Window's method for 
determining the distribution of electric field gradient in 
certain well known ferrous complexes, e.g. PeSO^.TH^O, 
Fe(]ffi^) 2 (S 0 ^) 2 . 6 H 20 , Na 2 [E’e(CN) 5 No] . 2 H 2 O (SKP) and an amorphous 
system [34] 35 H'a 20 + 55 Si 02 + 10 Fe 20 ^. The isomer shift 
distribution of K^Fe( 01 T)g. 3 H 2 O and 310 stainless steel have also 
been analyzed by this method* The present work introduces the 
viewpoint that this distribution method need not be restricted 
to amorphous systems only, but can be utilized successfully 
for crystalline systems, improving the agreement with the 
experimental data as compared to standard Lorentzian fitting 
and giving information about the local defects or environmental 
changes around the Fe-ion, 


3 . 3.2 Computational method 


For an BFG distribution p(lvl) we can in general 
describe the Q^. doublet line shape f(2) by a distribution 
p(lYj) of the Lorentzian lines in the form 


f(Z) 


h 


/ 


^ 03 


pdvi ) 

1 + 4[{Z>V)/¥]2 


d7 


(5.1) 


where V is the principal ociaponent of BP6, ¥ is the full 
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Lorentzian 'width at half maximum and h is the amplitude of the 
line. 

p(lvl) is expanded in a cosine series 'with the boundary 
conditions 


(-E) 


v=o 


(^) 


= 0 


v=v_ 


max 


p(lv|) = 


! [Oos {fSi) . (.i)*!] 
n=l max 


(3.2) 


Y is the cutoff value of V, taicen far away from the 

Iu3,X 

peak so that 




0 


The Pourier coefficients b^ are calculated by the best 
fit to the data. If IT -»■ ” any p(lvl) can be represented by 
equation (3.2). Since the coefficients b^ are calculated by 
solving (IT+I) simultaneous equations, the value of N is 
restricted by the amount of computer time one can afford and 
the efficiency of the computer, A value of W of the order of 
15 is adequate. However the authors have tried upto a value 
of IT = 40 in the DEO-10 Gomputer. 

For the isomer shift distribution the same method was 
used and hence all the above relationships hold good with Y 
replaced by isomer shift S. 
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5.5.5 Results and discussion 

The results of the BP6 distribution analysis for the 
spectrum of the well known standard substance sodium rdtro- 
prusside (National Bureau of Standards momber 725), obtained 
from Radio Chemical Center, Amersham, U.K. , ar« shown in Pigs. 
5.7 and 5.8. Pigure 5.7(a) shows the quadrupole splitting data 
compared with the best fit by a pair of lorentzian lines of 
arbitrary position, width and height. Pigure 5.7(b) shows 
the same data with best fit to equation (5.1). If can be seen 
that the experimental points are always slightly shifted 
outwords with respect to the Lorentzian curve where as file 
curve corresponding fo fhe EPG distribution (Pig. 5.7b) 
passes through the experimental points at the peak. The 
distribution curve p(lV(), wherefisp'ut in velocity units, is 
shown in Fig. 5.8, This fit is obtained with N = 25 and ^ e 
Q V = 6 mm/sec. The mean quadrupole splitting at the 
maximiim value of p(M) is 1.705 mm/sec which agrees well with 
the reported values [29,50], The PWH14 of the distribxition 
around the mean value is 0.12 mm/sec. The knee in the distri- 
bution curve is due to the presence of dehydrated sodium 
nitropresside present in the material. The position of the 
maximum of this minor peak is 1.970 mm/sec which agrees well 
with the quadrupole splitting 1.964 mm/sec of the dehydrated 
sodium nitroprusside [51]. The oscillations and the negative 
value of p( {v| ) af low and high values of V are because of the 
truncation of the Fourier seriesk 




VELCCITY Cmm/stcl 

'Saebaiier spectra of SHP (a) pure Iioreatsian leart squar 
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The EPG- distribution in Pe30^.7H20 is shown in Pig. 3.9. 
The mean value of ^ e QV as obtained from the peal position 
of the distribution is 3.25 mm/sec and the P¥HM of the distri- 
bution is 0.24 mm/sec. The distribution for Pe(lE[^ ) 2 (SO^) 2 . 
6 H 2 O is shown in Pig. 3.10. The mean value of the quadrupole 
splitting ^ e Q V is 1.732 mm/sec and the PWByi is 0.235 mm/sec. 

Mbssbauer spectra of the amorphous [35 Na 20 + 55 Si 02 + 
10 Pe 20 ^] system with the fitting curve, using equation ( 3 . 1 )? 
is shown in Pig, 3.111. The quadrupole splitting ^^sing 
lorentzian fit is 0.962 mm/sec. On careful examination of the 
spectrum it can be seen that the resultant spectrum is because 
of three overlapping quadrupole split doublets. The correspon- 
ding quadrupole splittings are precisely obtained from the 
distribution curve (Pig. 3.12), the values being 1.01 mm/sec, 
0.45 mm/ sec and 1.66 mm/sec. The first one is due to the 
Pe-atom at the tetrahedral site and the second one due to 
that in the octahedral site. The original of the third one 
is still unidentified and requires other confirmatory experi- 
ments. In a normal lorentzian fit the peaks were difficult 
to resolve beyond a doublet whereas the p(lvl) distribution 
gave clearly the three separate doublets. The octahedral site 
peak is the most dominant one, with PWHM of 0.41 mm/sec. 

Table 3.1 gives the consolidated information on the 
quadrupole splittings obtained from the lorentzian fittings 
as well as from EPG distribution curves for the samples 



P(V) ARBITRARY ‘UNIT 



Hi* 5*8 . If® ia f*Cs^4l 7H2O 

at »#©• 
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TABLE 3.1 


The quadrupole splitting obtained by lorentzian fit, and from 
p(|vj) distribution and PW'HM of the P(IVI) distribution. 


Substance 

Q.S, from 
lorentzian 
fit 

(mm/sec) 

Mean Q.S, 
from peak 
position 
of P(1V») 
distribu- 
tion 
(mm/ sec ) 

S tandard 
value in 
literature 

(mm/sec ) 

F.VJ.H.M. of 

P( ]V 1) dis- 
tribution 

(mm/sec) 

3KB 


1.704 

1.705 

1.7034 [29 ] 

^ 0.12 

Dehydrated 

SEP 

- 

1.970 

1.964[31] 

- 

PeS04.7H20 


3.200 

3.250 

3.217[32] 

0.24 

Ferrous 

Ammonium 

Sulphate (FAS) 

1.714 

1.732 

1.738[32] 

0.235 

So. ium silicate 
glass 





(a) Tetrahedral 

0^962 

l.OlO 

1.00 [33] 

0.41 

(b) Octahedral 

- 

0.450 

0.5[33] 

- 

(c) Unidentified 

— 

1.660 

— 

- 


^Numbers in square brackets indicates the references. 
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studied by us. It is evident that the distribution method 
gives more detailed information as well as better fit than the 
Lorentzian fitting method. The line widths in each case is a 
measure of the inhomogeneity of the EPG- across the volume of 
the sample. 

The isomer shift distribution of potassium ferrocyanide 
(K^I'e(CIl)g.3H20) is obtained using the above method and is 
shown in Pig. 3-14. The isomer shift as determined from the 
least square fit with a single lorentzian line (Pig, 3.15) is 
0.220 mm/sec. with respect to sodium nitroprusside and that 
from the position of the peak of the distribution curve 
(Pig. 3.14) is 0.21 mm/sec. The PWHM is 0,32 mm/sec. The 
isomer shift distribution of 310 stainless steel is shown in 
Pig. 3.15. The peak position is 0.175 mm/seo and the PWHM is 
0.42 mm/sec. 

The method presented here open up the possibility of 
looking into the structural origin of the SPG and isomer shift 
again and develop Mosaic models for deviations from the ideal 
structure in the hitherto accepted perfect crystalline 
systems. 


3.4 Recoilless Praction 

The recoilless fraction of y-^rajs in the Ferrous 
Complexes FeSO^. 7H20f Pe(NH^)2.(S0^)2.6H20, K^Pe(ClJ)g,3H20 
(PPG) and Ma2(Pe(GU)^lfa).2H20 were measured and the Debye 
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Isomer snift distribution in 310 stainless 
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temperature is calculated from this. Computer programmes 
were made for calculating these parameter from the Mossbauer 
data, using the following method. 

For a lattice in which all nuclei occupy equivalent 
lattice sites the recoilless fraction (Debye Waller factor) is 
given as 


f = exp(- <x^> 


(3.3) 


2 P 2 

where k is reduced wavelength of y-^ay ( = 0.0188 A 

5 7 2 

for 14.4 keV y-ray of Fe )5 <x is the mean squai’e displace- 
ment of the Mbssbauer atom in the direction of y-rsiy at tempera- 
ture T, averaged over the life time of the nuclear excited 
state (t = 0.97 x 10”"^ s for ^"^Fe). The Mbssbauer effect can 
be detected in systems satisfying the condition <x^>^ <<>: 
from Bqn. ( 3 . 3 ). 

Expanding <x in normal co-ordinates and using the 
Debye approximation for the phonon spectrum 

I/Sn 

f = exp[- ^ (1 + 4(^)2 / 2 E^] ( 3 . 4 ) 

^D 0 e -1 

2 2 —2 

where E is the recoil energy (R = B^/ 2 Hc = 2 x 10 “" eV for 
5 7 

Fe) k is Boltzman constant and 0 ^ the Debye temperature. 

The Debye Waller factor is calculated using the method 
in Ref. 2 by finding the integrated intensity (axea) of the 
Mbssbauer spectrum. 

The total area of the Mbssbauer absorption spectra 


is given by 
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t -V2 

ho-t = S’rtfgTle [I^(t/2) + l3_(V2)] (3.5) 

where I^(t/2) is hyperbolic Bessel function of order n and 

r 

T = (N f cr(0) ~— ) is the effective thickness of the absorber. 

Bl, i 

2 

N is the number of atoms/cm of the absorber and f recoilless 
fraction of the absorber, r is the absorber line width corrected 
for absorber thickness using the expression given by 0 'Connor 
[12] and Shirley and Kaplan [11] 

r = 2 r (1 + 0.135 t) for 0 < t < 5 or 

T = 2 r (1 + 0.145 t - 0.0025 t^) for 4 < t < 10 

is the natural line width of the y-^ay cj(0) is the resonance 
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cross section for Pe, The correction of experimental data 
due to non-resonant radiation was taken into account by the 
factor S. S was determined by the filter technique given by 
Honsely et al. [35]. A brass filter of thickness 0.005 in. is 
used with the same geometry to find out the background radiation. 
For a source with no high enex’gy Z-rays a filter of 0.005 inch 
of brass will allow less than 0.02 ’/, of the 14,4 keV gamma rays 
through, while removing only 4 /. of the high energy radiation. 

The count rate obtained after putting the filter between the 
absorber and detector multiplied by 1.04 is the desired back- 
ground rate due to high energy gammas. If is the 

number of 14.4 gamma rays and the background as determined, 


then 
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Knowing the recoilless fraction, the Debye temperature can 
be calculated using equation (3,4). In the limit of low or 
high temperatures compared to Eqn. 3.4 can be approximated 
to 


and 


f = 


f = 




6Ej^© 

exp (- 5 


) 


for © « ©jj 
for © » ©jj 


(3.8) 


Using the above described method the recoilless fraction 
of the various iron complexes ahe determined and is shovm in 
Table 3.2. 


The recoilless fraction is calculated using • the method 
of O’Connor et. al. [37,38] also as follows. The area under 
the absorption dip is given by 


A 




(3.9) 


where l(v) is the intensity observed when the source- 
absorber relative velocity is v. 


I(v) = I( “) (l-fg) + I( » ) fg. 

OJ t o'(Vt-v) IVt 

/ S(V 3 _) exp [- SfoT — ~ ^ (3.10) 

«« oa 

m 

where f is the recoilless fraction of the source, S(v) is 
s 

the normalised source energy distribution, t is the effective 
absorber thickness at its peak value and cf(v) is the nuclear 
absorption cross section 
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TABLE 3.2 


The recoilless fraction, Debye temperature and thickness 


corrected 

EWHM of some 

Ferrous Complexes. 


Substance 

Recoilless 

Fraction 

Debye 

Temperature 

K 

Thickness Corrected 
Width 
(mm/sec . ) 

SEE 

0.469 

252 

0.250 

EAS 

0.511 

246 

0.248 

PEC 

0.562 

200 

0.258 

EeS0^.7H2 

0 0.671 

519 

0.210 
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t = N f a(0) 
a 

Substituting these values 

CO CO t C5 ’(Vt ) 

A = / dv / S(v+Vj^) 1-exp [ — ] dv^ 

— CO — 00 

“ t cr(v) 

= f / l-exp [ — — ] dY 

® .« a(0) 

Expressing v in terms of B, the energy of the y-^ays 

A = fs / l-exp [- — dE (3.11) 

«. CO 

o(B) is assumed to have a Voigt profile as 

a(E) = a(0) — / 2 

(S-y)^+l 

where y = 2(B-E^)/r where r is the natural line width and P 
is a Graussian broadening parameter (equal to “ for pure 
Lorentzian). The area A was computed for different values 
of P and t and the measured areas are fitted against t ciirves 
to determine the values of t and hence the recoilless fraction 

f . 

a 

The area A of the absorption dip was measured by the 
following method which is similar to the ’black absorber’ 
method [35] for determining the recoilless fraction. In this 
case the black absorber is replaced by a ’white source’, which 
has a uniform recoilless intensity over a velocity range + V 
and zero intensity elsewhere. Then the fraction of radiation 
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absorbed by an absorber is directly proportional, to tlie integrated 
area, of the absorption dip within + V. The wbJite source is 
simulated by moving an actual source to scan the range + Y and 
integrating the total intensity transmitted in this range. If 
the source line width is negligibly small then the effective 
intensity per unit velocity interval, in the range + V when in 
motion would be proportional to (2V)“^ and the total transmitted 
intensity would be 

I(v) = I( - ) (1 - ^) (3.12) 


where A is the absorption dip area within the range + V, In 
actual source the deviation from the ideal 'white source' is 
less than 2 /. at 0.7 V for V equal to ten times the natural 
line width. As V increases the approximation improves and 
the proportion of the absorption dip area effectively spanned 
by the source approaches its saturation value. Measurements 
were made for Y between 40 to 300 natu- al line widths over 
which the graph of l(v) against (2Y)“^ was found to be linear. 
A was then found from this graph using the relation 


A 


iro 

d(2Y)"^ 


and I( ” ) was found by extrapolation of the straight line 
to (2V)”^ = 0. The values of f determined using this method 

Q. 

closely agrees with the value obtained by the method described 


earlier 
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CHAPTER 4 


ELECTRIC FIELD QRADIEHT AT IRON PURITIES IN hop METALS 


4.1 Introductiori. 


The nuclear quadrupole interaction in metals has been 
the subject of considerable interest in recent years. The 
origin of the electric field gradients (SPG) at the nuclei in 
metals is still not well understood. However, the extensive 
experimental data on BPG in non-cubic metals reported [1,2] 
in recent years has enabled us to draw some conclusions 
regarding the origin of BPG's, in paJ'ticular the relative role 
of the lattice contribution and electronic contributions 

(eq^l). In this chapter, we present a brief review of the 
observed systematics, and the theoretical and semi- empirical 
models for SPG determination and our study of the EPG at the 
iron impurity in some hep metals namely Ti, Se and Fd. 


4.1.1 The EPG tensor 


The electric field gradient tensor at the origin 
due to a point charge Ze is given by 



ft 



(4.1) 


where r^ = I x.^ and x- are the Cartesian co-ordinates of the 
i ^ 

charge. Being a traceless tensor only five of the nine compo- 
nents of ¥. - are independent. A unique axis system called the 
'principal axes of the BPC tensor* can he defined such that 
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the off-diagonal elements vanish and the diagonal elements ahe 
chosen so as 



The EPG tensor is usually described by only two independent 
parameters, = eq and the asymmetry parameter r) = 

- '^yy^'^'^zz* above ordering p is restricted to 

0 < r) < 1 . 

The EPG is mainly contributed by two sources [3]j 
(i) fixed positive ions of the crystal surrounding the 
Mbssbauer atom in non-cubic symmetry (the lattice contribu- 
tion), and (ii) the non-cubic electron distribution in 
partially filled valence orbitals of the Mbssbauer atom (the 
electronic contributioxi) , The resultant electric field 
gradient can be written as 

^<1 = * ^'Sei (4.2) 


where (l-Yn-) is the Sternheimer antishielding factor (for 
iron = - 9.14) a well defined value which 

can be calculated by the lattice sum method. In tho case 
of lattice with hexagonal close packed (hep) structure the 
EPG tensor is axially symmetric with the main axis along the 
c-direction and equal for all lattice sites. The value of 
qiatt is calculated numerically by Eas and Pomerantz [4] 
which is given by 


®%att 


[0.0065 - 4.3584 (| - (1)^)] (4.3) 

4'ae^a^ 
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The e<lg2 contribution is defined much less precisely and is 
difficult to calculate. Any evaluation of would require 

a detailed knowledge of the wave functions of conduction 
electrons in the impurity- host system. One would expect that 
®lel strongly depend on the anisotropy of the .Fermi 

surface, i.e. on the orbital character of conduction electrons. 
The spatial distribution of conduction electrons near the 
impurity should follow the symmetrj?' of the host matrix. This 
approach has been justifiably used by Watson et al. [3] in 
their calculation of eqQ]_ for pure metals and ordered alloys. 

4.2 Systematic Trends of the EPg’s in Metals 


4.2.1 Correlation between and 

Raghavan et al. [5] after studying the systematic 
trends of EPS in non-cubic metals pointed out that there 
exists an almost linear relationship between eq^^ ®9latt 
(1 -yP as 


ecigi = - K: (1 - tP (4.. 

where K:fr2-5 for moderate values of ®9q^-tt 

Por the cases where both the magnitude and sign of the 
BPG- are known experimentally, determined by subtrac- 

ting the theoretically calculated values of 
from the experimental data. ®9@i then plotted against 
®^latt (Raghavan’s universal correlation). However 

Raghavan et al. did not coaasider the rare-earth metals since 



9 




the unfilled 4f shell quite delinked with the lattice, is the 
dominant source of the electronic part of the EFS. Qd has 
been included because the configuration of ion is 

spherical, making 4f contribution to EFG equal to zero. On 
the other hand cases such as Te and Tel are omitted since 
the EPG in these cases is dominated by local covalency 
effects [6], 

The recent data on non-cubic metals are plotted in a 
similar way by Kaufmann and Vianden [1] and are reproduced 
in Fig. 4.1 with minor modifications. It is seen that most 
of these data support the correlation suggested between the 
electronic and ionic contributions. 

Recently, Krusch and Forker [7] have shown that there 
are cases in which either the sign or the magnitude of the 
electronic SFG are incompatible with this correlation. In 
all those cases studied the host metal o'” the probe atom 
belong to a transition element series which suggests that in 
transition element systems the electronic EFG may be dominated 
by contributions, which do not arise in simple sp-metals. In 
the case of scandium host and iron impurity the constant 
K 18 [7]. In a theoretical study Piecuch and Janot [8,9] have 
shown that the magnitude of the local contribution to the 
BFG at transition element impurities should be proportional 
to the occupation of the d-states of the impurity, which 
depends on the dens ity of states of the host. 



.g* 4»1* ^ eleetrimle imd lG33t»^ eoiEitrl 

bcitiG^ t&^-&L9''$M''hm94. m »9&mM t«tsa* 

Eigare tak«a fr®« 3P«f*[l3. 
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4 . 2.2 Dependence of EFG on the probe Used 

The experimental data on SPl in non-cubic metals shows 
that it is probe-dependent in the form of v.alence difference 
between impurity and host. Over and above the host electronic 
contribution an additional local electronic contribution may' 
be present from the valence electrons of the probe axom which 
remain spatially correlated to the host atom. An expression 
which reveals the possible influence of the impurity on the 
resulting BPG is the ratio of this BFG at an impurity site 
P in a host H to the BEG of the pure metal (P = H), i.e. 

(H,H) (4.5) 

It is evident that the BEG at an impurity site should be 
very sensitive to changes in the conduction electron wave 
functions around this site. These changes can be brought 
out by the electron redistribution around impurities by 
which the charge difference between impurity and host matrix 
is screened out within short distance around the impurity 
atom. To show that this charge screening manifests in the 
BEG, a plot of the set of ratios of Eq. (4.5) for the 
different probes P in each matrix H as a function of the 
valence difference between impurity and host 62 is repro- 
duced [ 10 ] in Eig. 4.2. Somewhat similar plot of Collins 
[ 11 ] is also reproduced in Eig. 4.2, which brings out the 
fact that in sp metals, the electronic BEG is approximately 
proportional to the valence of the impurity which in turn 
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is related to the number of conduction electrons screening 
the impurity. 

4.2.3 Dependence of the EPG on temperature 

The quadrupole interaction frequency (QIP) is found 
to vary with temperature. For pure metals and dilute 
impurity alloys the variation approximately obeys a simple 
empirical relation by Christiansen et al. [12] 

eq(T) = eq(0) [l - (4.6) 

where T is the temperature* and 3 is the slope parameter. 
Besides this expression another test function 

eq(T) = eq(0) [l-B^T-B 2 T^] , (4.7) 

where eq(0), B^^ and B 2 being free fit parameters, was also 
found suitable for many systems. 

The slope parameter B in Eq, (4.6) is positive 
excepting TaBe [13] and the magnitude of B is related to 
the strength of the coupling between the probe and the host. 
The values of B were found to be different for different 
impurities in the same host for some systems. Also, for a 
given impurity the values are different in different host 
matrices. In the case of rare-earth systems with partially 
filled 4f-shells the above relation does not hold good. 

Since most of the EPG contribution is likely to come from 
the 4f-shells delinked from the lattice, Kishiyama et al# 
[14] used a pseudopotential approach including the influence 
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of lattice vibrations to find the temperature dependence of 
tne EFG in non-cubic metals. The resulting eq(T) factorizes 
into a Debye-Waller factor and a lattice sum over screened 
ions as 

eq(T) = ^^“^eff^ ®^ion^^^ I (4.8) 

where is the lattice sum over screened ions, kp the 

2 

Permi vector and <u is the mean-square atomic displacements 
at temperature T. 

An exact explanation for the observed changes in the 
slope parameter B is difficult to give, since a large number 
of mechanisms are responsible for the observed temperatiAre 
dependence. While interpreting the data, one should exercise 
caution to observe that the changes are not because of some 
changes in the phase of the matrix, the lattice locatioxa of 
the probe atoms etc, 

4 . 3 Theoretical Estimates of BPG*s 


4.3.1 The lattice contribution 


The lattice contribution to the field gradient eq i: 


eq 


■latt 


= I 2e2 [ 


N 


5 Oos^e^f - 1 




■] 


where and Ojj being the polar position coordinates of the 

•4- l-n 

N ion with respect to the EPS principal axes and suimaation has 
to be carried out over all ions in the lattice except one at 
the origin. Eijboer and de -Wette [15] used the Fourier 
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transform metiiod to calculate the lattice sum. This -sum has 
been evaluated numerically by T.P. Das et al. [4] to get the 
value for an hep lattice (Bq. (4.3)). 

4.3.2 Electronic contribution 


'The electronic contribution to the EP& at nuclear 
site in metals can be written as 

= - e / [-'" — -3- — ] (i't 

T ^ 

where is itie conduction electron density, the calcu- 

lation of which is complicated due to the lack of exact 
wave functions [16] p(r) is given by 

p(r) = 2 

nk ■ 


where the summation over n refers to bands and that over k 
to the states in eaoh band that are occupied, that is those 
lying below the Fermi surface. The factor of 2 takes care 
of the two spin states for eaoh electron. Using the above 
expression for p(r) 


^el 


= -e 


nk 


"'Pnk 


(?) 


Oos^e - 


1 


o 


nk 


(?) 


(4.9) 


The Wave functions tpjjg(?) could be obtained using appropriate 
crystal potential and OPW method. 

Recently Thompson et al. [17] Rave presented a first 
principle analysis of the contributions from the isotropic, 
anisotropic and aJaharmonic lattice vibrations to ifee lattice 
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and electronic SFG. The temperature-dependent pseudopotential 
was incorporated to account for the temperature dependence. 

It is found that the lattice SPG is primarily affected by the 
anisotropic phonons, but not by the isotropic vibrations. 
However, the isotropic component almost entirely accounts for 
the variation of the electronic SPG with temperature. 

4.3.3 Pseudopotential approach and charge screening ef. tacts 

The pseudopotential method can avoid the difficulties 
in the calculation of the exact crystal potentials, by 
constructing explicit electron wave functions which may be 
used to evaluate the EPG’s. The basic method in these 
calculation is to find the electi’onic charge distortion 
which causes a new charge distribution. An impurity ion in 
a metal having an impurity-host valence difference results 
in a, redistribution of the conduction e) "jctrons around the 
impurity to screen the effective charge differences. 

Nishiyama et al. [14,18] have suggested a basically 
different model for PPG. This model is based on the simul- 
taneous treatment of electronic screening along with the 
lattice vibrations. The effect of conduction electrons is 
to shield the valence ions so that their effective charge 
becomes 

Ze 6(r-R) - e (r-f) (4,10) 

p is the screening charge distribution, 
s c 
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The EFG- at the origin due to a single screened ion at 
posi'tion. R with respect to the probe is given by 

F(R) = / [Ze 6(r-R) - e D^^(r-S)]r^ 2?. (OosO^)d^r 

The net EFG- at the probe nucleus is evaluated by 
taking contributions from all the screened ions antishield od 
by the screened probe ion. A new enhai^cement factor 
is also used in these calculations. So the time dependent 
EFG at temperature T is given as 

eq(t,T) = (1-Yeff)/^' 6(E-R^(t)) F(R) d^R 

i 

^ th 

in which Rj_(t) is the instantaneous position of the i ion. 

12 

Since the frequency of the lattice vibrations ( 10 Hz) 

O 

is large compared to the QIF ( < 10 Hz), we have to take 
the time average over the ionic vibrations which are regarded 
approximately as uncorrolated and isotro'^ic [19], so that we 
get the result 

eq(T) = eq(t,T) = (1-Yeff)/ T 6(S-H. ( t )) F(R) d^R (4.11) 

The Fourier transform of this equation leads to 

eq(T) = / (S^(k)-l) exp[-(?.x)^/3] F(k) d% 

.... ( 4 . 12 ) 

Where F(k) is the Fourier transform of F(R) and S(k) the 
structure factor of the non vibrating lattice. The influence 

of the lattice vibrations on the BFG at the probe ion is 

2 

described by the Debye— kaller factor with <z > being the 
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isotropic mean square displacement. 3y the use of a simple 
pseudopotential V(R) derived from a Coulomb potential V*^(R) 
screened by a free- electron gas 

V(fi) = -n f 2^ g . jrP n 1,2 
2n2 i if H 

equation (4.12) may be finally written as 

eq(T) = (1-Yeff) 3 4 (4.13) 

where = P ”^ 2 (^ 1 ) 2 P 2 (CosOj^) over the screened ions, 

and kp is the Permi radius. e(k) is the dielectric function 

[ 20 ]. 

This model has not given results consistent with 
experimental data for the magnitude of SPG- but the tempera- 
ture dependence agrees with that of many systems. 

4 . 3.4 Conduction electron charge shift model of Bodenstedt 
and Perscheid 

Bodenstedt and Perscheid [21] suggested a model for 
the origin of EPG in hep metals essentially in terms of the 
charge shift of the conduction electronsj this shift is 
related to the deviation of from the ideal value of 

3 . 

^ = (-2)2 for an hep lattice. They started from the assump- 
tion that the conduction electrons fill the space between 
the ions. For an easy analysis the conduction electron 
distribution is separated into charge clouds, the centres 
of which are situated in the middle between nearest 
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neig'h.’bour ions as shown in Pig, 4.3. Each positive ion is 
thus surrounded hy six negative charge cloiids in its own 
hexagonal plane and three clouds above and below. The latter 
ones form planes of pure negative charge in the middle 
between subsequent hexagonal planes of ions. The chaJ'ge 
clouds in hexagonal plane are on equivalent sites and must 
therefore bear the same charge Qj^p. The same ar,gum 9 nt holds 
for the charge clouds between the hexagonal planes, let that 
charge be called Since the tot.oL lattice is neutral 

%p ^ '^el = " 5 ^ 

where © is ibe effective charge of the ion. The conduc- 

tion electron charge shift 6 is defined by 

% = " i ^eff ® 

%1 = " ^ ^eff ® 

By using the approximation that each of the conduction 
electron charge clouds and the ionic charge distribution 
have spherical symmetry, a lattice sum calculation can be 
done for each of them yielding 

, eq = eq^^^ - I " I (4.14) 

In order to estimate the charge shift 6 for a given 
metal lattice with a given q/a, one can start from, measured 
elastic coefficients and calculate the stress needed in 
order to deform the id^al h<^ lattice into the shape with 



t c-axis 



conduction electron charge clouds 


Fig. 4-3 Charge distribution in on hep metal idealized for the 
Chorge Shift model. Figure token from ref. C213 
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observed c/a. Then one can calculate the charge shift £ which 
produces this stress by the electrostatic forces. For a 
hexagonal close packed lattice the value of 6 obtained in 
this way is given by 


c/i 


(I)' 




^2 2 
^eff ® 


3' 


(4.15) 


(-23^^ + S 


12 


Sll) 


+ S 


33 


13 


] 


are the elastic coefficients of the host lattice. 
For comparison with experimental data one has to 
apply the Sternheimer correction factor (I-Yq,) as 


since all charges contributing to eq are outside the probe 
ion sphere. 

The effect of the vibrations of the probe nucleus 
is incorporated by calculating the average value of eq('i^),'r 
representing the instantaneous position of the probe ion. 
Then one gets 

2 

eq(r) = eq(0) + l/2l — <x| >(l-e) 

3X, 

2 

<xi > 

with e = — ^ 

<4 > P 

2 3 eq. 

The quality ^ = I — 2^ lattice sum method 

3x| i 

I 24 P.(Cos9. ) 

^ 1 V 

for the positive ions and conduction electron charge clouds. 
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This model of conduction electron charge clouds also 
allow us to estimate the probe dependence of the for probe 
ions with an ionic radius not larger than that of the lattice 
ion. To incorporate the probe effect, first the total lattice 
sum in pure metal eq is determined and the EFG due to the 


twelve surrounding charge clouds carrying the effective charge 

P 

of the host ion ,eq is subtracted. Then the EPG of the 


twelve surrounding charge clouds with the effective chaJ'ge of 
P 

the probe eq (Z' is added. In order to get the total SFG 


this is m'oltiplied with probe. 


e 


Total 

'^Probe 


(l-'Yc^)[eq-eq^(Zejf) + (4.l6] 


is calculated using the lattice sum method as 


P 1 ^eff ® , r,, 1 (-1 I (c/a)^ - 3 


'0 a- 


[i + i (c/a)"]- 


The results of such calculations for zinc are reproduced in 


Fig. 4.4. 

The different models which are discussed have only 
limited success in predicting the EPG in metals. Experimental 
determination of EFG' s in different systems will be valuable 
in studying the various systematic trends which enables us to 
have a better understanding of the origin of the resultant 
EFG observed in metals. However , o'Ur experimental results 
and discussions wo-uld be better uiiderstood in view of these 


discussions. 
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Fig 4-4 The dependence of EFG in Zinc on the probe 

used. The solid line represents the prediction 
of the Chorge Shift mode!. Figure token from 
ref. LZn 
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tedru pole Inter ations at in TltariiT:m Metal 

4.4.1 Introduction 

The main objective -of the work reported here is to 
measure the electric field gradient experienced by ^'^Fe 
probe doped in Titanium metal as a dilute impurity, using 
i^bssbauer spectroscopy. The results are discussed in the 
framework of the available theoretical and semi- empirical 
models. The distribution of EPG at the probe nuclei is also 
determined from which am estimate of the different locations 
in which the iron atom has gone is found out. 

Titanium is a transition metal in 3d-series with a 
hexagonal close packed structure. Piecuch and '^anot [9] 
have calculated the BFG- and its temperature dependence in 
transition metals of the 3d, 4d and 5d series using a first 
order perturbation dynamical theory in random phase appro- 
ximation. The anisotropy of the metal v/as treated as a wea,k 
perturbation to the electronic Hamiltonian with respect to 
a cubic environment. The EFG obtained in this way in pure 
titanium is (- 7.6to+2.5) 10^^ V/cm^, 

Three measurements of the quadimpole splitting in 
^"^PeTi at room temperature are available in literature, 

Quim [22] reported a fairly large value for the quadrupole 
splitting pBp = 0.32 + 0,03 mm/see, while Wortmann and 
Williamson [24] got a very sflSall splitting of 0,1 + C,02 
mm/sec through their Mhssbauer experim^ts. Both of them 
us8d the method of electroplatii^g Co acti¥dt;y on Ti foiXs 
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4.4 Quadru^ole Interations at in Titarii'um Metal 

4.4.1 Introduction 

Ihe main objective -of the work reported, here is to 
measure the electric field gradient experienced by ^'^Fe 
probe doped in Titanium metal as a dilute impurity, using 
tfossbauer spectroscopy. The results are discussed in the 
framework of the available theoretical and semi- empirical 
models. The distribution of EPd at the probe nuclei is also 
determined from which an estimate of the different locations 
in which the iron atom has gone is found out. 

Titanium is a transition metal in 5d-series with a 
hexagonal close packed structure. Piecuch and ‘^anot [9] 
have calculated the EFG and its temperature dependence in 
transition metals of the 3d, 4d and 5d series using a first 
order perturbation dynamical theory in random phase appro- 
ximation. Ihe anisotropy of the metal was treated as a weak 
perturbation to the electronic Hamiltonian with respect to 
a cubic environment. The EPG obtained in this way in pure 
titanium is (- 7.6 to +2.5) 10^^ Y/cm^. 

Three measurements of the quadi^pole splitting in 
57 

EeTi at room temperature are available in literature. 

Quim [22] reported a fairly large value for the qu^rupole 

splitting aEq = 0.32 + 0.05 mm/ sec, while Wortmann and 

Williamson [24] got a very small splitting of 0.1 + C.02 

mm/ sec through their Mbssbauer experiments. Both of them 

. 57 

used the method of electroplating Oo activity on Ti foils 
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and diffusing at 1200°C. Devere and De ¥aard implanted the 
Co activity in a single crystal of Ti at an energy of 
130 keV and got an upper limit of = 0.12 mm/sec by 

I'ldssbauer measurements, after annealing at 700^0 for a period 
of 1 hour and assuming that the ^"^Go activity has gone to 
unique substitutional sites. 

Since the above experiments with source matrix as 

titanium ,gave such widely varying results we have done the 

measurements, by preparing an absorber of pure titaniiim with 

57-m 

Ee as a very dilute impurity for the first time, and using 
57 

the Co in Rh matrix as the source, 

4.4.2 Sample preparation 

I'itanium of 99.99 /. purity and iron ('^0,2 at 7* ) 

57 

enriched in Pe isotope, were melted in a fused alumina 
crusible in an induction furnace at high vacuum. The melting 
Was repeated a couple of times to get a hemogeneous alloy. 

The alloy thus obtained was powdered with a diamond file and 
Mbssbauer measurements were done with this unannealed sample. 
The quadrupole splitting of this sample was measured at room 
temperature. The resultant spectrum obtained showed a quadru- 
pole doublet spectrum of appreciable line broadening. 

The alloy was then annealed at 800°0 for 6 days after 
sealing it in a quartz tube under high vacuum and slowly 
cooled to room temperature. This alloy was powdered and the 
samples were prepared by keeping the powder in a copper ring 
with cdlophane tape with a distribution of ^0 mg/cm of powder 
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and also by mixing the powder with the aTaldite adhesive 
followed by drying in the form of a tliin plate. This method 
ensux es the homogeneity of the sample and helps the sample 
to be perfectly random in orientation. 

X-ray diffraction patterns of titanium and of the 
prepared alloys were recorded and compared. The analysis 
snowed that the alloy has the same crystal structure as that 
of pure titanium, 

4.4.3 Experimental 

The details of the experimental set-up were given in 
Chapter 2. The data were collected for sufficiently long 
time to ensure good absorption peals and the effects of 
the statistical fluctuations were small. Pfeiffer et al. 

[25] and Zobayashi et al . [26] pointed out that if the peaks 
are not well separated with overlapping broad lines and the 
statistics is poor, an independent fit of Lorentzian para- 
meters to the data may lead to erroneous results. A para- 
meter £ = Dyz Was set for statistical quality of the spoctmm 
where D is the amplitude and N is the number of counts per 
channel at' the peak. It was shown that independent fit 
procedure is valid only if 5 >10. In all our measurements 
this parameter was kept well above this limit. 

4.4.4 Results and Discussion 

(a) Quadrupole splitting at room temp^ature 

The Mossbauer paTaiieters obtained in our measurements 
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is shown in Table 4,1. The average Q.S. at room temperature 
[Pig. 4.5a] in the unannealed sample is 0.343(8} mm/sec. An 
electric field gradient of 1,63x10^'^ V/cm^ is deduced assuming 
the asymmetry parameter to be zero. 

ihe annealed sample M'dssbauer spectra [Pig, 4.5b] 
clearly showed the presence of two different well defined 
sites for the iron atom corresponding to Q, splitting of 
0.112(8) mm/ sec or to an EPG of 0,53(4 )xlO^'^ V/cm^ for the 
inner lines and 0.380(8) mm/sec corresponding to an EPG of 
1,86(4 )xl0^'^ V/cm^ for the outer lines. 

Pollowing the arguments of Devare et al. if we assume 
that the smaller quadrupole splitting corresponds to the 
substitutional site of the iron atom the value of Q.S, 
obtained agrees quite well with the value reported by 
Wortmann and Williamson, Moreover this method of sample 
preparation adopted reduces very much the chances of iron 
forming clusters. Comparing the resiuLts with the reported 
interme tallies of Ti with Pe [27] it is most probable that 
the other site of iron in our alloy is interstitial which 
gives a quadrupole splitting of 0.39(1) mm/sec. The agree- 
ment of the observed line widths (0.275 mm/sec) with the 

57 

expected values further shows that the Pe has gone to 
unique substitutional and interstitial sites in Titanium 
matrix for the absorbers prepared by us. Purther, after 
preparing several samples it is found that getting a unique 
substitutional site alone for ^'Fe in the host matrix 
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TABLE 4.1 

Observed Mbssbauer parameters for ^^Pe^ system. 


Temperature Quadrupole isomer Shift* eq(l+p^/3 ^ 

°K ^^m/seS^ mm/sec in 10^'^ Y/om^ 


80 

0.128 

(6) 

0.172 

(7) 

0.613 

(28) 

125 

0.126 

(6) 

0.160 

(8) 

0.602 

(28) 

295 

0.112 

(8) 

-0.016 

(5) 

■0.538 

(58) 

560 

0.105 

(6) 

-0.16 

(8) 

0.491 

(28) 

586 

0.098 

(8) 

-0.21 

(9) 

0.468 

t 

(38) 

'With respect to iron metal. 

Numbers 

in ; 

paranthesis 

show 


the errors. 



COUNTS 
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titanium is very difficult with the usual methods of alloy 
preparations. 

b) Comparison with EPG-' s a,t other probes in titanium 

Since titanium is an hep metal the field gradients 
are expected to be axially symmetric. Therefore we have 
analysed our data assuming p = 0. Using the known value of 
the quadrupole moment of the 14.4 keV state of ^'^Ee (0.192 b) 
and the measured q^uadrupole splitting, the electric field 
gradient eq is obtained. The ionic contribution is calculated 
using the Das and Pomerantz [4] formula 

^'llatt = [0.0065 - 4.5584 (| - 1.635)] (4.17) 

and the known Sternheimer antishielding connection factor 
(1 -Yq(j) = 10.14 for Pe. The ionic contribution for 

^'^EeTi system is 0.47x10^"^ V/cm^. If we assume that the signs 
of and eq^^ opposite as observed for most of the 

other systems, the electronic contribution eq^^^ is obtained 
to be 1.01x10^^ V/cm^. 

Table 4.2 shows the EFG at different probes in titanium 

57 

host reported so far. It seems that Pe in titanium produces 

the minimum EPG at the impurity site if we assume that the 

minimiim Q.S. observed corresponds to a substitutional site. 

On the other hand if the larger Q. splitting is taken as the 

17 

one corresponding to substitutional site the IPG is 1,86x10 
V/cm^ which is comparable with the BPG obtained with other 
probes. Since the atomic voltime of Pe is comparitiTeiy roaller 
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(7.1 compai’ed to 10.6 for Ti) there is enough possibility of Fe 

going int erstitially in a diffusion process. Considering the 

Hume— Rot her y rules for solid solution formation -A— is to be 

AK 

less than 0.15 and in this case for Fe and Ti -v 0.15 which 
is not very much lower than 0.15, The above studies indicates 
that further ion channeling experiments have to be done to 
locate the exact site of the Fe atoms before arriving at a 
conclusion regarding the substitutional position of iron atoms 
in titanium host. The calculated value [(-7.6 to 2.5)xlO^^ V/cm^] 
of EFG in pure titanium by Pieuch and Janot [9] is much smaller 
than the measured field gradients in titanium. Table 4.2 also 
shows the electronic contributions to the ERG obtained by 


substacting the ionic contribution from the total EFG, The 

quantities eq^^^ and eq^q are the electronic parts obtained by 

taking the signs of and eq^q to be same and opposite 

57 

respectively. The sign of EFG at Fe probe in titanium is 
negative [23]. The sign of EFG at Od and HF impurities are 


given as positive [30,32]. But the ^q^q calculations reveal 
that if the resultant* BFG is positive at Gd the ratio of 
®9el/®'^ionl ~ which deviates much from the Raghavan's 

systematics [5] where as if it is assumed to be negative the 
ratiojeq^q/eq^Q^t^omes to 2.14 which agrees well with the 
systematics . 


c ) Correlation between sqqQjj ®^@1 

Table 4.3 gives the measured BFG*s in the hop transition 
VQ.lL JL3.I3-L 0 

metal hosts/ so far in literature. The electronic contributions 
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to the EFG- in hep transition metal hosts are plotted against 
the ionic contributions in Pig, 4»6, The quantity K = 
has a value of 2,15* This shows that the astern ^"^PeTi agrees 
well with the systematic correlation between and 

Considering the other site corresponding to an observed EPG- 
of 1*86x10 V/cm if the sign of this is positive^ the ratio 
of eq^^^^ I comes to 2, 95, and if the sign is negative 

d ) T emp er at\ir e d ep end enc e 

From the measurements done between room temperature 
and liquid nitrogen temperature it is seen than this system 
almost obeys the law for the variation of SPG (Pig. 4.7) 

as 

eq(T) = eq(0) [1 - B 

The parameters obtained from these measurements are 
eq(0) = 0.63 X 10^ V/cm^ and 

B = 32 X 

The experimental findings of Q.S, are compared with the 
predictions of the charge shift model, later in this Chapter 
(Table 4.7). 

e) EPG distribution at the impuriiy site in titanium host 

Using the method described in Chapter 3 the distribu- 

57 

tion of EPG at the two different sites of Fe in titanium is 




'Fjf.4'7 Temperature dependence ot EFG in ^'^PeT 
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found out, and is shown in Pig, 4.8. The width of the distri- 
bution for the site corresponding to the smaller value of Q.S. 
is 0.012 mm/sec and that for other site 0.016 mm/sec. 


f) Goniparison of experimental data with conduction electron 
charge shift model 


The charge shift model calculations described earlier 
are done to find out the EFG in titanium and its temperature 
dependenc e. 

The compliance constants are determined from the known 
values of the elastic constants 0. . of titanium at different 


temperatures as 


S 


11 


^11^33 " ^13 


(G11-C12) (0^^833 4 - C 3_2033 - 20^3) 

.2 


G 


S 


12 


13 


^12^33 


^^ll-^lT^ ^^11^33 ^12^33 ’2^13^ 


-C 


13 


’13 
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^ 11^33 ^ 12^33 
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33 
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Table 4.4 gives the values of a and c/a of titanium at various 

temperatures. The values of ®%p ®^el taken 

from Ref. 21 for the known values of c/a for titanium. Table 

^eff ^ 

4.5 gives these values in tmits 5 

4iEe^a^ 

nature or c/a value* 


as a function of t^pe- 
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TABLE 4.4 


Lattice parameters of titanium at different temperatures. 


Temperature 

a 

i 

c/a 

4 

2.9416 

1.5843 

73 

2.9435 

1.5850 

123 

2.9449 

1.5856 

293 

2.9496 

1.5874 

386 

2.9526 

1.5878 

523 

2.9558 

1.5882 

656 

2.9593 

1.5886 

774 

2.9634 

1.5894 

920 

2.9670 

1.5941 


Data taken from Ref. [36] 
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TABLE 4.5 


The value of ^'^hp ®^el different c/a values of 

titanium at different temperatures in units of qIA%z^s? . 


Temperature 

Og. 

c/ a 

^lion 

®^hp 

^'Jel 

4 

1.5843 

0.2195 

-39.566 

41.701 

73 

1.5850 

0.2164 

-39.597 

41.691 

123 

1.5856 

0.2158 

-39.600 

41.682 

293 

1.5874 

0.2059 

-39.647 

41.656 

386 

1.5878 

0.2041 

-39.658 

41.650 

523 

1.5882 

0.2024 

-39.668 

41.644 

656 

1.5886 

0.2006 

-39.679 

41.638 

774 

1.5894 

0.1971 

-39.699 

41.626 

920 

1.5941 

0.1765 

-39.823 

41.556 



125 


The Values of the elastic constants of titanium at 

different temperatures are shown in Table 4.6. Table 4.7 shows 

the calculated values of EFG at ^'^Fe probe in titanium host 

as a function of temperature for different values of Z'^.. 

eft 

The Sternheimer anti shi elding factor (l-v ) is taken to be 

cc 

57 

10.14 for the Pe probe in these calculations* 

In the above calculations we have neglected the effects 
of the vibrations of the probe ions. ‘ Unlike the Uishiyama et al. 

[18] model, the isotropic vibrations of probe ions do not 
contribute to the EFG in charge shift model. The contributions 
from the anisotropic vibrations may be expressed as 

2 

eq( vibrations) = ^ (l-e) <X^| > 

2 2 2 2 

where s = <Xj_ > /<X^^ > and Xj_ and Xj| are the mean squared 
displacements of the probe ions perpendicular to the c-axis 
and parallel to the c-axis respectively 

There are two factors which mainly contribute to the 
temperature dependence of EFG. The contributions come from the 
variations of the elastic constants and the lattice parameters 
with temperature and from the anisotropic vibrations of the 
probe ions. The temperature dependence of EFG is quite sensi- 
tive to the anisotropic vibrations of the probe ions, 

57 

■ The measured values of the field gradient at ‘Pe in 

titanium host agrees well with the experimental values for 
2 

= 1 aiid Z* ^ = 2* The charge shift medial calculations 
^ eff ^ 

i ed to other hep Bctals and the results are shown along 
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TABLE 4.6 


The values of elastic constants of titanium at different 

12 ? 

temperatures in units of 10 dynes/cm . 


Temperature 

H 

H 

O 

Gi 2 

^13 

G 35 

4 

1.761 

0.869 

0.683 

1.905 

73 

1.749 

0.871 

0.68 

1.894 

123 

1.726 

0.877 

0.681 

1.876 

298 

1.624 

0.92 

0.69 

1.807 

373 

1.579 

0.934 

0.694 

1.774 

523 

1.495 

0.961 

0.692 

1.715 

673 

1.416 

0.978 

0.69 

1.661 

773 

1.368 

0.985 

0.688 

1.627 

923 

1.299 

0.992 

0.688 

1.567 


Data taken from Ref. [57] 
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TABLE 4»7 


The values of electric field gradients at ^"^Ee in titanium host 
at different temperatures calculated using the charge shift 


model. The 

value of is 

taken as 4. 



Temperature 

(°K) 



^eff 


.... - D /"s'v* 4 

0 

1 

1.5 

2 


4 

0.375 

0.502 

0 . 566 

0.63 

- 

73 

0.364 

0.484 

0.552 

0.615 

0.61 

123 

0.361 

0.483 

0.545 

0 . 606 

0.602 

293 

0.335 

0.448 

0.504 

0.561 

0.538 

386 

0.325 

0.435 

0.489 

0.544 

0.468 

523 

0.313 

0.418 

0.470 

0.523 

- 

656 

0.300 

0.401 

0.451 

0.502 

- 

774 

0.290 

0.387 

0.435 

0,484 

- 

920 

0.253 

0.337 

0.3787 

0.420 

- 
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with the experimental values in Table 4.8. 

4.4.5 Conclusions 

The electric field gradient at doped substitutionally 

in the transition metal titanium as a dilute impurity is 
measured to be 0. 538(38 )xlO^'^ V/cm^ at room temperature. The 
results of the measurements of EPC for ^"^PeTi agrees with the 
charge shift model calculations. The ratio of ®Q. 0 i/®l 3 _on 
agrees well with the universal correlation relation of Raghavan 
et al. 

*57 

4.5 Quadrupole Interactions at Fe in ' Trigonal Selenium 
4.5.1 Introduction 

In the last few years considerable progress has been 
made in the understanding of the physical properties of selenium. 
It is known that selenium exists in the amorphous, hexagonal 
^ ,trigonal) and a- and p-monoclinic forms. The growth morpho- 
logy and physical properties are reported in detail in references 
[38,39], Trigonal selenium with its chain lattices and strong 
anisotropy have a number of peculiar and interesting physical 
properties. Since selenium is a saaiconductor, it will be 
interesting to study the BPG produced at an impurity of the 
transition metal group, and to find out the electronic contri- 
bution to the BPS, because of the more localised nature of the 
conduction electrons. 
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4.5.2 Crystal structure of hexagonal selenium 

A molecule of hexagonal selenium consists of atomic 
chains displaced in Zig-Zag fashion in the direction [0001] 

(Fig. 4.9). The shortest interatomic distance in the chain and 
the shortest distance between the chains of hexagonal selenium 
are 2,32 E and 3.53 £ respectively. The covalent bond angle 
Se-Se-Se is 103'^. Between the atoms in the rings and chains, 
in the case of hexagonal selenium, there exist homopolar links 
and the rings and chains themselves are bound to one another 
by considerably non-polar forces. Hexagonal selenium belongs 
to the space group 'S'5^1 with unit cell dimensions a = 4.3662 A 
and c = 4.9536 A with p = 120°. The unit cell has three 
molecules lying on a helical chain of symmetry 3^ ox 32* These 
chains pass through the Corners of a hexagonal lattice. The 
atoms lie on special positions on a two-fold rotation axis 
so that only one position parameter x, jeded to be determined. 

4.5.3 Sample preparation 

trigonal selenium crystals were grown by the method 

described by G.H. Griffiths and B. Fitton [40], Selenium 99.99 7. 

57 

purity and iron powder (--^0.1 at */. ) enriched to 91 A. in 'Fe 
were sealed in a quartz tube under high vacuum. The sealed 
samples were kept at 600°C for about 10 days, after which the 
ampoules were transferred to a furnace kept at gN0°O, for 4 days, 
in which case it crystallizes in the form of spherulites. This 
Was cooled to room temperature and the ^herulites were powdered. 
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The X-ray diffraction pattern was obtained using the XRD-5 X-ray 
diffractometer to determine the crystal structure. The results 
of X— ray analysis shows that selenium has grown into the trigonal 
form with the unit cell diamensions given earlier. About 30 mg 
of powder was sandwitched between two cellotapes covering an 
area of nearly one square cm. The thickness of the active 
material was therefore about 0,05 mg/cm , which was quite small 
so that the broadening due to sample thickness could be neglected. 
A separate sample was prepared by thoroughly mixing the powder 
specimen with araldite adhesive for getting a distribution of 
pure random orientation. This sample was used for the measure- 
ments from liquid nitrogen temperature to room temperature. 

4,5.4 Results and discussions 

Since the melting point of selenium is 217°C, the range 
of temperatures at which the Mossbauer spectra can be studied, 
is limited to comparatively low temperatures. 

^'or pure electric quadrupole interaction, the quadrupole 
splitting in the case of may be written as 

2 p. „2 1 

= 4^ (1 + f 

¥liere eq is the electric field gradient at the nucleaX* site and 
T] is the asynmetry paraisieter* Values of eq and t) cannot be 
determined separately from Mbssbauar spaotra alone and 

therefore one has to obtain the value of rj from some other 
experiments to obtain the matgnxtude of idio electric field 

- n ^ ISKcs'jai’hjali Al* 
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with I obtained a value of r| = 0.5 + 0.02 for trigonal 
selenium. 

A typical M'dssbauer spectra obtained by us is shown in 
Fig. 4.10. The Mdssbauer parameters derived from the data are 
given in Table 4.9. The values of isomer shift are given with 
respect to iron metal. Using the known quadiupole moment 
Q = 0,192 barn of the 14.4 keV state of ^"^Fe, the experimental 
values of eq(l + p /3)^ were calculated. 

Earlier experiments [41,42] by Boolchand et al. in 
trigonal selenium yielded a value 110x10^ V/cm^ for the BFG 
with ^^^Te in Se, and 69 . 2x10^ V/cm^ with . In these 

125 

experiments the line widths are -^6 mm/sec, in the case of Te 

and 2 mm/sec in the case of I. In the present experiment 

from the observed line widths and* solubility data, it appears 
57 

that Fe is substitutional in selenium matrix. The BFG 
obtained with the above values for p ana Q at room temperature 
is eq = 0. 749 ( 9 )xl0^'^ V/cm. Boolchand et al. [41] reported recen- 
tly that the sign of BFG on in selenium is negative. At 

present it appears to be difficult to comment on the sign of 
BFG for ^"^Fe in selenium. The lattice contribution to the BFG, 
®%on ~ 0-763 ^ 10^'^ ?/cm^, using Baa and Pomerantz relation 

[Eq. 4.17] and with (i-Y^) ioT ^'^Fe equal to 10.14. Using this 

^ €0 

value we obtain = 1.512, if the sign of BFG is positive 

and sign of eq. negative, and eq ^ = 0.014 if the sign of EUl 
is positive and sign of also positive. That is, in 

+V. contribution becomes negligible and 
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TABLE 4.9 

Observed Mdssbauer parameters of ^'^F eSe . 


Tempe- 

rature 

K 

Quadrupole 

Splitting 

(mm/ sec) 

Isomer 

Shifts 

(mm/sec) 

TT 

eq(l + ^ )2 
in 10^'^ V/cm^ 

eq 

in 10^'^V/cm^ 

82 

0..175 (6) 

0.450 (5) 

0.837 

(28) 

0.804 (27) 

123 

0.170 (6) 

0.451 (6) 

0.813 

(28) 

0.781 (27) 

195 

0.167 (8) 

0.456 (7) 

0.799 

(38) 

0.767 (56) 

298 

0.163 (7) 

0.467 (6) 

0.779 

(33) 

0.749 (32) 


%ith respect to iron metal. 
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all the El’S is of ionic origin, fable 4.10 gives the BEG for 
different probes in selenium host. It j.s seen that besides the 
ionic part, the values are highly probe dependent. This 

is not unexpected because, the co- valent bonding between the 
atoms on the same chain play an important role in determining 
the EFG, and the bonding of the probe atom with neai’est host 
atoms could change the EEG considerably. 

a) Correlation between eq^j and 

The present data on selenium seem to qualitatively 
support the Raghavan' s correlation if we assume that and 

eqei are of opposite signs. 

b) Temperature variation of QIE 

The temperature dependence of the QIF's [14,17,43] may 
be expected from the changes in lattice parameters, from the 
lattice vibrations, or possibly from th.. interactions of the 
neighbouring quadrupoles. Since the number of observations 
between room temperature and liquid nitrogen temperature is 
small, it is difficult to establish the actual relation for 
the temperature dependence of QIF. However, we checked that 
our data fit with empirical relation 

eq(T) = eq(0) (1 - B 

c) Recoilless fraction 

57 

The recoilless fraction of FeSe absorber is measured 
by using the method adopted by Trooster et al* [44] with some 
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TABLE 4.10 

The currently available data on quadrupole interactions for 

I II 

different probes in Selenium host, are 

calculated values assuming the total BEGr to be positive and 

17 2 

negative respectively. All the values of EEG are in 10 V/cm . 


Probe 

®^exp 

*11011 

I 

®1el 

, II 
®1el 


+110.0®- 

1.45 

108.55 

-111.45 

129j 

+ 69.2 

1.34 

67.86 

- 


+ 0.749^ 

0.763 

- 0.014 

- 1.512 


(a) Reference 1 (b) Present work. 
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modifications. The method described in Ohapter 3 can be applied 

only if we know the number of active nuclei /cm)of the 

to determine 

absorber which is difficult/in this case. So the following 
method is used to determine the recoilless fraction (f ), 

fg^ is related to the mean square displacement of the 
Mbssbauer atom in the absorber by the expression 

= exp <x^>^] (4.18) 


and ilnf is proportional to the absolute temperature T for 

oL 

T/O^ >0.5 as 


- m f (T) 

9, 


6RT^ 


(4,19) 


where R is the recoil energy of the Mbssbauer atom? the 
Debye temperature and k is the Boltzman constant. In other 
words 

- £n f (T) = -XT 

f (T) is also given by the equation 

3. 

^ 

o 

where t is the effective thickness of the absorber, N the 

3 

p 

number of Mbssbauer atoms/ cm and 0^ is the resonance absorp- 
tion cross section for the Mbssbaner radiation [Chapter 3], 


an t = £n Ho + in f (T) 

Q, 0 3 

5= an Ho* + X T 

P 
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at temperature 


tn = in + V Tj_ 


and at temperature T2 


Xn t = an No„ + X 

3-2 02 


from which 


an = 


To t - Tt an t 


T2 - T^ 


( 4 * 20 ) 


t and t are calculated from the Mossbauer spectra areac at tv 
^X ^2 

temperatures. The area of an absorption peak can be calculated 


from the transmission integral as 

2 f 


ir(v^) = B [i-B ri- 


(V^)' 


S ^ S^ /’tt'ttN 2 , / p / ^ \ 2 


Ti^o (¥-Vg)^ + {'^^/ 2 y 


X exp {- t 


(r /2)' 


^ + (r /2)^ 


>dV]] 


(4.21) 


where N(V ) is the number of y-quanta that reach the detector 

B 


when source and absorber have a reJ.ative velocity with 

s 


respect to each other, P is given by 

O 


^s 


^■0 = ■ 

® K + Kioo) 


which is the recoileas fraction of the incident radiation 
corrected for background. H^Cos) is the number of y-quanta 
reaching the detector pa“ unit time when the relative velocity 
between source and absorber is very larg^» is 'the back^ound 
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rate when a 0.005 inch brass plate is introduced between the 
absorber and detector for the same geometry. F can also be 
determined from the Mossbauer spectrum of a substance like 
SKF or natural iron foil where the distribution and the 
recoilless fraction are known, is the natural line width 
of the source and Y is the energy variable expressed in equi- 
valent Doppler velocity. 

For small values of t the exponential can bo approxi- 
mated by 

N(¥ ) = B [1 - F . ■■ 5-] (4.22) 

s ® ^ V + (r/2)^ 

s 

where r = 2 . The peak intensity and absorption area are 

given by 

= 1/2 F^t and = 1/2 ti F t 

o s a o ' s a o 

Trooster et al. [44] made a table of t^ for various values of 

cl 

area A, of a Lorentzian line fitted to N(V_) as given in Iq. 

s 

(4.21). These values were obtained by calculating N(Y_) 

O 

numerically and subsequently fitting a Lorentzian line to the 
calculated spectra. If the source line width deviates from 
by an amount as a first approximation the measured line 

width should be corrected by this amount, before it is used for 
the determination of t , i.e. (r - 4r„) should be used instead 

9 p 

of r . 

57 

Using the above method the recoilless fraotion of 'Fe 
in Se is determined at room temperature and the avwage value 
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comes out to be = 0,18. The Debye temperature of the probe 
atom in the host lattice can be determined using Eq. 4.19. The 
mean square displacement of the host is related to that of the 
probe through the relation 


2 

<x > . 


m 


hostij- 


^ for T 


<x > 


host 


■21 

imp 


0°K 


where is the mass of the host atom and is that of 

the impurity atom. The Debye temperature of the host can be 

2 


calc’ulated from <x > host using eq. 4.19. 

The Debye temperature of the probe atom in the host can 
also be determined from the slope of, logarithm of the effective 
thickness against temperature graph. Prom Eq. 4.19 the slope 
of this straight-line graph is 


6 E. 


£n t . - an t 


a- 


3,1 


kG 


D 


Ti - T2 


from which 9p can be fourd out. Using the above method the 

57 

calculated value of Debye temperature of Pe in selenium is 

120°K. 


The author has also applied this method to find the 

57 

recoilless fraction of Pe in the isoelectronio host tellurium 
and the value of f thus obtained at room temperature is 0.43. 

3 


4,5.5 Conclusion 

The electric field gradient at doped in the seailcon-- 

ductor selenium as a dilute impurity is measured to toe 0,749(43)2 
17 2 

10 ¥/ciii . The electrojaic eoatrilmtioii to the total BPS is 
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negligible if the experimental and ionic EPG's are jtaken as 
positive. Otherwise it obeys the universal correlation of 
Raghavan et al. 

4.6 Quadrupole Interation at ^'^Ee in the hop Rare-Earth 
Reodymium as a Dilute Impurity 

4.6,1 Introduction 

Reodymium is a rare earth metal of the group Illb 
series with a double hexagonal close packed structure, 
which means that the layer sequency is ABCABC instead of ABAB, 
Therefore two inequivalent lattice sites with differing 

by 10 /. should exist [45]. Raghavan et al. [5] showed that 
the EFG in rare earth metals does not obey the universal 
correlation proposed by them because the electronic part of 
EFG is mainly contributed by the unfill d 4f shell. Recently 
Ernst et al. [45] proposed a new classification for hexagonal 
metals. They have shown with the available data at present that 
the electronic contribution to the total EFG is proportional to 
the ionic part, the proportionality constant K being positive 
for all group Illb and IVb hexagonal metals and negative for 
all group lib, Vllb and Vlllb metals. In rare earths the 
studies are mainly concentrated on the heavy rare earth metals 
from Gd to Lu with more than half filled 4f shell where one 
expects the conduction bands to be quite similar. These hop 
metals have lattice structures with almost the same o/a ratio. 

To check the validity of these proposed systematics for a light 
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rare-earth we have taken the system ^"^PeNd and measured the SPG 
by Mbssbauer effect. 

4.6.2 Sample preparation 

Neodymium (99.9 1 . purity and iron (--^0.2 at T', ) enriched 
57 

in Pe were melted in an induction furnace in Zirconia Crusible 
under high vacuum. The resulting alloy was sealed in a quartz 
tube under high vacuum and annealed at 500^^0 for 4 days. Since 
neodymium oxidises very easily in atmosphere the alloy, after 
breaking the seal, was kept in mineral oil. This shining alloy 
ball Was powdered with a dimond file in mineral oil and later 
mixed with epoxy resin to prevent oxidation and to have a unifoni 
distribution. The quadrupole splitting of this sample was 
measured at room temperature and at liquid nitrogen temperature. 

X-ray diffraction pattern of this sample was studied and 
the analysis showed that this alloy has the same crystal struc- 
ture as that of pure neodymium and no other intermetallic 
compounds are formed in significant proportions. 

4.6.3 Results and discussions 

57 

Ihe typical Mbssbauer spectra of the FeM system at 
room temperature is shown in Pig. 4.11. The observed Mbssbauer 
parameters at room temperature and liquid nitrogen temperature 
are shown in Table 4.11. The quadrupole splitting at room 
temperature is 0.30(1) mm/ sec which is an average value obtained 
from many independent runs. An electric field gradient of 
1.45(5 )xlO^'^ Y/cm^ is deduced from this, assuming the aspmaetjy 
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TABL3 4.11 

c: n 

Observed Mbssbauer parameters of B eM . 


Temperature 

Og 

Quadrupole 

Splitting 

mm/ sec 

Isomer Shift* 

mm/sec 

eq(l+Ti^/3)^^^ 
in 10^^ V/cm^ 

295 

0.304 (9) 

-0.140 (8) 

1.454 (43) 

175 

0.317 (8) 

0.06 (7) 

1.516 (38) 

82 

0.331 (8) 

0.209 (6) 

1.584 (38) 


*V/ith respect to iron metal. 
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parameter to be zero. The PWEM of the 'absorption lines observed 

0.32 mm/sec which is comparable with the natural line width at 

57 

room temperature of Co source in rhodium matrix and natural 

57 

iron absorber, which shows that the Fe has gone to unique 
positions, most probably substitutional, in the alloy. The 
isomer shift at room temperature obtained is -0.140(8) mm/sec 
with respect to iron metal and it changes to 0.209(6) mm/sec 
at 82 K. 

The lattice contribution to the EFG by using the Das 
and Pomerantz relation for hep metals is 0,08x10 V/cm which 
is very small compared to the experimental value. This shows 
that the electronic contribution is the major source of BPG in 
this case and the ratio of is 1 ^ 17 . The systema- 

tics proposed by the Ernst et al. [45] » gives for heavy rare 
earth metals, a value for ! ©'lel^'^ionl ^ value of 
and © 9^2 with this systematic behaviour is reproduced in Table 
4.12. Our results of F eM shows that the electronic contri- 
bution is very much augmented for the less than half filled 4f 
shell . 

4.6.4 Some theoretical aspects of EFG in rare earth metals 

Pelzl [46] has given some insight into the role played 
by the conduction electrons and the 4f electrons in heavy rare 
earth metals. The total EFG is because of three sources. 1. 

The EFG produced by the 4f electrons, 2, The EFG arising from 
the point charges surrounding the rare-earth ions, and 3, The 
EFG produced by the conduction electrons. Bach of these 
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TABLE 4.12 


Systematica of electric field gradientsin pure hexagonal group 
Ilh, Illb (including the rare earths), IVb, Yllb, and Vlllb 
metals. The experimental EFG data the iorAc contri- 
bution = ®^latt^^“’^co^ ^ electronic part eq^^ = 

“ ®q-i are listed. All the values are given in units 

XOXi 

of lO-^ ' 7/cm^ . Predicted signs for eqgjj.p> whose signs have 
not been measured, and the corresponding values for eq^^ are 
put in parentheses. 


Group 

lattice 

®^exp 

®%on 


Ilb 

Zn 

+3.0(2) 

-1.96 

+5.0 


Gd 

+6.7(9) 

-3.65 

+10.5 


Gd (24 kbar) 

+6.0(9) 

-2.90 

+8.9 


Gd (50 kbar) 

+5.4(8) 

-2.60 

+8.0 


Gd (74 kbar) 

+4.8(8) 

-2.50 

+7,3 

Illb 

Sc 

(+)0.38(2) 

+0.22 

(+0.16) 


La 

(+)1.5(2) 

+c .50 

(+1.0) 


Gd 

+2.8(3) 

+1.06 

+1.7 


Tb 

+3. 4a. 7) 

+1.32 

+2.1 


Dy 

+4. 3(1. 5) 

+1.53 

+2.8 


Tm 

+5.0(3) 

+1 * 64 

+3.4 


Lu 

+3.6(5) 

+1.36 

+2.2 

IVb 

Ti 

( + )1.2(1) 

+0.40 

(+0.8) 


Hf 

+9.5(4) 

+2.84 

+6.7 

Vllb 

Tc 

(-)0.7(4) 

+1.56 

(-2.3) 


Re 

-4.9(2) 

+2,08 

-7.0 

Vlllb 

Pe (e phase) 

(-)0.26(5) 

+0.44 

(-0.70) 


Go 

- 0 , 29 ( 2 ) 

+0.15 

-0.44 


Ru 

(-)0.49(14) 

+1.41 

(-1,90) 


Os 

-4.5(3) 

+3«24 

-7.7 


Bata taJcen fro® Ref, [45] 
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gradients are, modified aspherical distributions of the closed 
electron shells (Sternheimer effect). The resulting E?G is given 
as 

1 = (1 - Rq) + (1 - Rq ) q^jjCint) 

+ (1 - ^liatt 

gradient produced by the 4f electrons, l^att 
due to the positive ions in the neighbourhood, q^^Cint) and 
q^jiCext) are the contribution of the conduction electrons 
inside and outside the atomic sphere. The Sternheimer factor 
and account for the induced distortions 

The two field gradients can be related by introducing a 
parameter 6 as 

From a non-relativistic APW calculation ^as and Ray [48] 
obtained 6 = 0.4 for Dy metal. A nearly free electron behaviour 
is a good approximation for the conduction electrons , and a 
vanishing q^^Cint) is expected in heavy rare-earth metals. The 
electric field gradient produced by the 4f electrons which are 
well located inside the atomic sphere can be calculated using 
the relation [46] 

(1 - R^) q^^ == -3 e< J }ai J> 

< j| - 1/3 

the effective radial integral is defined as 
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<r-3> (1 . a. ) 

■where Rq is the Sternheimer factor which accounts for the 
contributions from the induced distortions of all other core 
shells due to the aspherical 4f electron distribution in the 
metal. An empirical relation is proposed which fits with the 
radial integral as <r^> = (Z - for 64 £Z < 70. 

For rare-earth metals the best fit yields = 1.940 and = 
58.41. 

From these theoretical consideration Pelzl [46] found 
that in the heavy rare-earth metals, the EPO due to lattice 
and conduction electrons should be due to the point charge 
contribution and an enhanced electron charge density at inter- 
stitial sites respectively, and is given by 

^ ^latt ^ ‘llatt 

It is also found that the conduction eleo crons also enhance 
the field gradient at the 4f site, 

4.6.5 Conclusions 

57 

The SPG at Pe as a dilute impurity in the light rare- 
earth metal neodymium was determined for the first time as ’ 
0.50(l)xl0^'^ V/cm^. The lattice contribution is 0.08x10^"^ 
V/cm and the ratio of jeq^^/eqj,-^^^|= 17, which violates the 
systematics proposed by Eaghavan et al. [5] and also' of Brnst 
et al. [45]. The IPG in this oase may be originating maialy 
from the less than half filled 4f shell. A. systematic study 
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of the EPG- in other hep rare-eafth metals with less than half 
filled 4f shell (Pr and Pm) will give a better insight to this 
behaviour. 


4.7 Oomparisons 

The three hep structures chosen to serve as the hosts 
57 

of Pe, represent three different catagories. Ti is a transi- 
tion metal of the 5d-series with c/a = 1.5873 and has interes- 
ting properties. The sign of field gradient in Ti depends on 
probe, For ^PeTi it is negative where as for “^OdTi it is 
positive. The ionic and electronic contributions follow the 
general trend observed by Raghavan et al. Selenium is a 
semiconductor which is isoelectronic in nature with Tellurium 
and Sulphur, The observed BPG has the same value as the ionic 
contribution to BPG (if eq is tahen as positive) which 
suggests that the conduction electron cc itribution is practi- 
cally zero in this case. Neodymium is from the rare-earth 
group with a less than half filled 4f shell and the results 
of our study shows that the 4f electron contribution is very 
much predominant over the other contributions to the total BPG, 
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CHAPTER 5 


HYPERFIHE lETTERAGTIOITS AMP PHASE TRAH3ITI0HS IN 
BIHARY IROH ALLOYS 

5.1 Introduction 

In recent years Mossbauer spectroscopj’’ has been found to be 
very useful in the investigation of order-disorder transitions 
in alloys and its effect on hyper fine magnetic field and 
quadrupole interaction. The broad distribution of magnetic 
hyperfine fields in the disordered alloys is notably reduced 
by atomic ordering, and the change in hyperfine field is explains 
as the result of the change in Po moment which increases with 
the increase in the number of Pe nearest neighbours. 

Of all the ME parameters, is the most sensitive 
for the microscopic changes in the structure of the alloys, and 
so the variation of H^^ with temperature and heat treatment, and 
its distribution because of the varying environment in the 
sample, will give a better understanding of the metallurgy of 
the alloy. 

Metal atoms coagulate some times in the non crystalline 
form which is in an extraordinarily disordered state i.e. in 
an amorphous state. The spectral lines of the hyperfine pattern 
of ferromagnetic amorphous alloys are very broad* At present 
there is general agreement that the broadening is oaoeed' by the 
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hyperfine field distribution* The poor resolution in the 
spectra led the scentist.s ’ to apply different models for 
the interpretation and evaluation of the spectra. There are 
in general three approaches [1,2,5] using the basic elements: 
1. microcrystallites, 2. molecules, molecular units, chanical 
compound clusters, embryos, 5. dense random packing. Kemeny 
et al. [4] suggested a ’’quasi-crystalline structure based 
on a locally distorted off-stoichiometric lattice'*. 

The presence of microcry staliites in these alloys 
seems unlikely because it would require that all the micro- 
crystallite arrangements exiiibit essentially the same spectra 
independent of the constituents. More over in contrast to 
experimental findings one would expect the spectra should be 
composed of a broad contribution, from the resonance atoms 
being located in rather distorted in between regions,and a 
rather sharp contribution of 1he atoms vj.^th in the grains. 

In the case of molecular models also one faces the 
same difficulty to envisage the solid as a whole with the 
molecules in their manifold of arrangements and orientations. 
The dense random packing model (Bernal model) seems suitable 
for fitting the spectra. 

The amorphous materials can be divided into two 
principal classes: 

i) the metglass alloys with o<Mapositlon iypically of om 
transition metal (Fe, Co, Ni) and one or several metal- 
loids (B, C, P, Si .). The composition usually 



alloy shows an interesting behaviour closely related -i-ath the 
order-disorder transition which appears at 600 As the 

long range order parameter of the fee based LI 2 type ordered 
structure of 1iie Pe^Pt alloy increases from 0 to 1, the 
martensitic transformation changes from the burst type to the 
thermoelastic type [7]. Since Pd lies between li and Pt in 
the VIII column of the periodic table the study of Pe-Pd alloy 
will be useful with respect to martensitic transformations. 

At high temperatures, iron and palladium form a 
continuous series of binary alloys with an fee lattice. At 
temperatures below 1100 K, however, a preference for surroun- 
ding with dissimilar neighbours becomes pertinent in the 
concentration range from 50 to 80 /. Pd. Around the PePd 
composition the alloy adopts the 11 q structure, and around 
the PePd- composition the II 2 structure. The crystal struc- 
tures as well as the regions of existence are shown in Pig. 

5.1a. The fcc-bcc martensitic transformation was reported 
[8] for an alloy of this type containing less than 20 at */, 

Pd. Recently Sohumura et al. [9] reported a thermoelastic 
fcc-fct martensitic transformation of Pe-Pd alloy which has 
not been observed before. The transformation is of the first 
order, but the axial ratio of c/a of the fet gradually decreases 
as the temperature lowers, vrhich is typical of the second 
order transition. The shape memory effect is also reported in 
this alloy. 

In order to elucidate the general features of the 

lattice transformation fcc-=i> 11 ==> fet, some details of the 

0 

olectron-miGroscopic observations in CuAu-X by Herabayashi 
et al. will be of interest [10]. Xn an early stage of the 
ordering treatment itself, smaU tetragonal crystallites of 
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the ordered phase come into existence. These ordered crystalli- 
tes are formed coherently with the disordered matrix, and due to 
the different lattice parameters an elastic strain energy 
develops. Since the c-axis contracts and the a- and b-axis 
expand during the fee ==^fct transformation, the elastic energy 
can be reduced when the ordered regions preferably form thin 
lamellae, lying in {110} planes of the fee matrix with the 
tetragonal c-axis along one of the cubic (100} directions. 
Nevertheless, some elastic stoain will build up around the 
tetragonal regions and this leads to micro twinning (Nig. 5.1b). 

In this process the c-axis of the tetragonal crystallites changes 
direction regularly between tvro of the throe cubic axes, thus 
reducing appreciably the strains in the plane of Nig. 5.1b, 
Perpendicular to this plane, however the strains remain, these 
are accomodated in a more advanced stage of the ordering by 
the interposition of the ordered material with the tetragonal 
c-axis in this direction. The atomic rearrangements with in the 
disordered phase preceding the transformation to the tetragonal 
structure remain undetermined. According to Zones [11], one 
might expect the tetragonal distortion to arise from the aniso- 
tropy of the ordering, and therefore a fair amount of ordering 
should occxir within the fee phase before the fee ==>fct trans- 
formation takes place. Hence a systematic study of the fee 
phase during an ordering treatment will be of coiasida’abl© 


interest 
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5.2.2 Sample preparation 

The alloy of Ee-^l.S at */, Pd was prepared by melting 
99.9 /. pure iron and 99.9 /. pure Pd in an induction furnace. 
This alloy was homogenised at 1100 for 5 days. This alloy 
X'/as powdered with a diamond file and annealed at 900 °G for 4 
hours and quenched into iced water. This powder was mixed with 
araldite adhesive resin so as to form a thin plate with a 

p 

distribution of 50 mg/ cm of Pe-Pd alloy, for M'dssbauer 
measurements, iUiother absorber was also made in a copper ring 
with a distribution of 50 mg/cm which is kept in position by 
celotapo for comparison. 

X-ray diffraction pattern of this alloy was recorded by 
the X-ray diffractrometer and the structure agrees xv'ell with 
that reported for the same alloy by Sohumura et al« [9] and 
Matsin et al. [21]. The diffraction peak corresponding to 220 
is shovm in Pig, 5.2. 

5.2.3 Results and discussions 
a) Results 

The observed Mbssbauer parameters for the Pe-31.2 at 
Pd alloy arc shown in Table 5.1. The hyperfine field distribu- 
tion were analysed by means of the procedure developed by 
Window [13]. In this, model independent evaluation, the observed 
line shapes are represented by a continuous pixsbability P(H) 
for the distribution of hyperfine fields, as 
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TABLE 3.1 


The observed Mbssbauer paraneters of Ee-51.2 7'. Pd alloy at 
different temperatures. 


Temperature 

K 

=ef£ 

EOe 

Q.S, 

(mm/sec ) 

Isomer Shift* 

(mm/sec) 

113 

307(2) 

0.44 (1) 

0.35 (1) 

126 

302(1) 

0.42 (1) 

0.31 (1) 

145 

299(3) 

0.32 (2) 

0.30 (2) 

169 

500(3) 

0.21 (2) 

0,26 (2) 

173 

299(2) 

0.19 (1) 

0.26 (2) 

198 

255(3) 

0.16 (2) 

0.23 (2) 

295 

299(2) 

0.11 (2) 

0.19 (1) 

389 

266(1) . 

0.11 (2) 

0.16 (2) 

430 

257(3) 

0.09 (2) 

0,16 (2) 

470 

206(1) 

0.09 (2) 

0.14 (2) 

508 

167(2) 

0.04 (2) 

0.06 (1) 

527 

99(1) 

- 

0.05 (2) 

540 

86(2) 

- 

0.06 (2) 

549 

39(1) 

- 

0.06 (2) 

600 

0 

- 

0.05 (2) 


%ith respect to iron metal* 
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P(H) = ! b„ [Cos SiS - {-1)=!] 

n=l ^ax 

with the boundaJ'y conditions 


( d P(H) v 

^ dH ^V=0 


^ dH 




0 . 


The hyperfine spectrum with this distribution can be then 
expressed as the modified Lorentzian form 


S(E) 


6 D P{H) dH 

2 j 

i=l l+4[(B-a^H)/¥^]2 


(5.1) 


where the spectrallines are assumed to be Lorentzian with 
position a^H, width W , and height D . . E is the resultant 
doppler energy corresponding to each point in the spectrum, 
is taken far away from the peaks so that 


The Fourier coefficients b^ are calculated by the best fit 
to the data. Since the coefficients b^ are calculated by 
solviiag a set of H+l simultaneous equations, M cannot be too 
large. With the computer programme developed by the author 
data was analysed with I of the order of 25, 

The Mbssbauer spectra of Fi^Pd alloy at various 
temperatures with fitting using Ihe distribution prograsaae is 
shown in Fig, 5,5. The corresponding hyperfiaa field distribu- 
tion is shown in Fig, 5.4, When the tssperature is sear the 
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Ourie "teiaperature the spectra appears to be of a broad single 
line which fits into a single line using the Lorentzian fitting 
programme. The data is again analysed using the isomer shift 
distribution programme described in Chapter 3 and the analysis 
shows clearly the hyperfine split spectrum at 549 The 

variation of H^^j(T)/H^^^( 295) with T is shown in ilg. 5.5. 

The paramagnetic Curio temperature determined from this graph 
is 554 The isomer shift in Fe-Pd is found to change 

with temperature. The variation of which is shown in Fig, 5.6. 

The temperature dependence of magnetisation of Fe-Pd 
was determined in residual field using the PAFG vibrating 
sample magnetometer. The curie temperature, T. of the sample 

w 1 

is determined using the interpolation programme. The value 

i 

obtained for T. in this way is 515 °K. The variation of i 

magnetisation a with temperature is shown in Fig, 5.7. The . 

value of T^ obtained agrees well with the extrapolated value | 

of Fe-Pd alloys studied by Matsui et al. [14], 

In order to study the fee ==> fet martensitic transfor- 
mation the M'dssbauer spectra of the sample was taJeen down upto 11' 
°K, The quadrupole splitting is found to increase from 0.07(2) 
mm/sec to 0.44(1) mm/sec as the temperature is decreased from 
room temperature to 113 The M^sshauer spectra of the sample 

shows that the fee phase is still present even ‘below the 
reported [9] transition temperature 173 asad at UJ it 
almost completely changes to the fet phase* ®he variation of 
Q.S, with temperature is shown in Table 5*2. fhs ooryaspoading 
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TABLE 5,2 

The quadrupole splitting of the fct phase of I’e-31.2 '/ , 
Pd alloy at different temperatures with the correspon- 
ding c/a values. 


Temperature 

K 

c/ a 

Q.S. 

(mm/ sec ) 

113 

0,965 

0.44(1) 

126 

0.97 

0.42(1) 

145 

0.975 

0.32(2) 

169 

0.99 

0.21(2) 

173 

0.995 

0.19(2) 



170 


c/a values from ref. [9] is also shown in the same table, 

B) Discussions 

The analysis of the M'bssbauer spectriom at vai^ious 
temperatures shows that the fee austenite phase of Be-Pd trans 
forms to a fet structure as can be seen from the changes in 
the Q.S. at around 173 However the spectrum also shows 

the coexistence of the fee and fet phases, at low temperature. 

(a) Temperature dependence of EPG 

The temperature dependence of the quadrupole inter- 
action with temperature [Table 5.2] is fitted vdth the empiric 
relation 

eq(T) = eq(0) [1 - S T^^^] 

A fitting of this equation with the experimental values of 
the fet phase gave a value of 

B = - 0,32 X 10“^ 

the value of eq(0) = 3.48 x 10^ Y/cm^. 

The dependence of BPG in Fe-Pd system on temperature is 
shown in Fig. 5.8 where the observed BFG’s are plotted 
against 
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(b) Universal correlation of EFG in tetragonal lattices 

The total EFG in tetragonal alloys is proportional to 
the lattice contribution to the EEG [15]. The temperature 
dependence of EEG in the fct phase of Ee-Pd shows that it 
varies almost linearly with the variation of the c/a ratio. A 
least squares fit of the function 

= h (=/a - op 

to the experimental points of the fct phase yields O 2 = 1*01. 
This agrees with the vanishing quadrupole splitting in the 
case of cubic symmetry. The variation of EEG for different 
values of c/ a is shown in Eig, 5.9. Eig, 5.10 sho^fs the 
variation of the lattice parameters and the c/a ratio for the 
Ee-31.2 t Eh alloy at various temperatures. 

(c) The distribution of hyperfine field and its temperature 
dependence 

The spectrum of Ee-Pd shows the following characteri- 
stics. Numbering the lines from negative to positive velocity, 
the following relations holds ^ Eg, T^> r^> 

effect is even easier to see in the different depths of the 
lines, this pattern can be understood from the model given by 
Window [16]. The distribution of the number of similar and 
dissimilar neighbours in the first shell causes a distribution 
of isomer shifts and hyperfine fields. The results of the 
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distribution studies is summarised by the values of H and AH, 
the position of the main peak in the P(H) distribution and its 
full width at half maximinm respectively, which is given in 
Table 5.3. It is seen that the most probable field coincides 
with the mean field within error limits in all cases. The 
width of distribution goes on increasing in the fee phase 
with increase of temperature which siiggests a more disordered 
phase. Whereas at low temperature side the width admost 
remain constant. 

In a disordered binary alloy, the Fe atom finds itself 
in different nearest-neighboxa- and next-neighbour environments. 
Hence in the disordered alloj' the Fe msirnetic moments have a 
distribution giving rise to a distribution of hyperfine fields, 
A complete random distribution is expected in the fee phase, 
since the order parameter is zero. At low temperatures, the 
distribution produces a spread in the saturation values of 
the hyperfine field [17]. If there exists a comiaon ordering 
temperature for all ions and if the temperature dependence of 
the hyperfine field for all ions is identical, then H should 
follow this temperature dependence and one expects that ^‘H 
shoxuLd steadily decrease as one goes from low temperature to I,,. 
Timoyoshi et al. [18] pointed out that the thermal -.verage of 
the magnetic moments on which the vreaker exchange field is 
acting decreases more rapidly x\'ith increase of temperatxxre 
than the one on which the stronger exchange fields. Therefore 
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The mean field most 

F¥HM (ah) of P(H) distri 
function of temperature. 


Temperature 

K 

^eff . 
kOe 

295 

299(1) 

389 

266(1) 

430 

237(3) 

470 

206(1) 

508 

167(3) 

527 

99 U) 


540 


36 (2) 


5.5 


probable 

field (H) and the 

ution in 

Fe-31.2 */, Pd as a 

H 

AH ^ 

kOe 

kOe 

298 

53 

257 

62 

216 

68 

197 

69 

142 

71 


112 


75 
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if th.e exchange field has local fluctuations at voxious lattice 
sites ; the reduced magnetisation has a broad distribution at 
high temperature. The range of the exchange interaction is a 
very important factor in considering the fluctuations of the 
exchange field. If the exchange interaction is a very long 
range one, the local fluctuation of the exchange field should 
be small even if the exchange interaction between different 
kinds of atom pairs has large differences, so that the magne- 
tisation and hence the internal, field have no distribution 
at high temperatures. On the other hand if the inteiaccicn 
is a very short range one, the -.l-loy should have a 1-rge local 
fluctuation of the exchange field according to the neighbouring 
atomic configuration. Thus the variations in the near 
neighbours environment also cause a distributiorx in exchange 
fields tbxTough out the material, and so, different i nes nay 
have hyperfine fields that show quite different teni..eratur'3 
dependences. This will have a tendency to cause -1 to 
increase as the temperature increases and H will lollov; some 
average temperature dependence. It is difficitit to cecide 
which of the above mechanisras is dominant. The appo?.reuce 
of two peaks in the p(H) distribution with two diffsrent 
values of H is an indication that there are two kinds of iron 
atoms with two different average magnetic moments end two 
different mean values for the hyperfine field at this temporatui'e. 
The appearence of the peak in the P(H) distribution centered 
H = 0 along with the peali with non-zero H value near the 


near 
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Gurie temperature tells us tha.t ferronagnetism does net vg,nish 
suddenly in these 8l.loys and also the ezistence of non-magnetic 
iron in the alloy [19]. Some ir“on atoms which are we?z3.y 
coupled by exchange interaction go to p-aramagnetic state at a 
comparatively lower temperature than the remaining iron atoms. 

The hyperfine field £.(0) at was found by extrapo- 
lating H to 0 °K and the valv-e obtained in this way is 514. k^e. 
The asjmmetry in the vjidths and heights of the peaks in the 
Mbssbauer spectrum of Fe-Pd alloy is due to the combined 
presence of the distributions of both q’uadrupole and dipole 
interactions. 

(d) Hecoilless fraction and Debye temperature of the Pe-Pd 
alloy 

The recoilless fraction is determined using the methc d 
of Trooster et al. described in Chapter 4. The averaps value 
of the recoilless fraction at room temperature ct the alloy 
is 0.588 and as the temperature is increased the value 
decreases to 0.42 at 470 It is interesting to note that 

for the fet phase the recoilloss fraction gees on decreasing v;ith 
decrease of temperature. The Deoye temperature at room tempera- 
ture is calculated from the r.jcoilless fraction as 279 ^4. 

The recoilloss fraction of the 44;- rl -alloy 
at various temperatures is slionmi in Table 5.4. 
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TABLE 5.4 

Recoilless fraction (f ) 3,nd effective t:hickness (t ) 

Si d. 

of the fee phase of Ee-51.2 y. Pd as a function of 
temperature. 


Temperature 

K 

t 

a 

f 

a 

295 

2.305 

0,558 

389 

1.975 

CD 

0 

0 

430 

1.792 

0.434 

470 

1.729 




5.2.4 Conclusion 
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The fcc==> fct martensitic transformation in the Pe-31.2 
at 7* Pd has been studied by the Mbssbauer analysis. The Q.S. 
clearly shows an incree.se after the trans forma cion. The 

variation of the EPG- in the fct phase obeys the universal 
correlation between lattice contribution and electronic contri- 
bution. The variation of Q.S, with temperature obeys the T^*^ 
law in the fct phase. The distribution of hyperfir.e fields 
has been studied by the model independent method of analysis. 
The width of the distribution increases with increase of 
temperature which reveals that the range of exchange inter- 
action is shorter in the Fe-Pcl alloy. The development of the 
paramagnetic peah even before the bulk curie temperature is 
reached shows that there is a distribution of curie temperatu- 
res in the Pe-Pd alloy studied. 


5 . 3 Mbs s bauer and I-iagnetic Studies of the AQor-jho\; s 




5,3.1 Introduction 

During the past twenty four years the class of 
compounds formed between rare-earths and 5d transition metals 
r.as gained particular interest, since one can beniiit l^ere 
from the intrinsic properties of both components, i.e. from 
the high magnetic moment per atom and the strong single-ion 
magnetocrystalline anisotropy of the rare-earth partner’, and 



181 


from the high magnetic coupling strength of the moments of the 
3d- transition metal partner. A subclass is formed by those 
compounds in which the rare-earth partner bears no magnetic 
moment (La, Lu, Y and titravaJ.ent Oe). Because of ti.e use in 
mal:ing permanent magnets much work has been done in the 
investigation of intermetallic compounds between rcre-earths 
and transition metals [20,21]. Strong interest in o-'orphous 
rare-earth (R)-l'e alloys began ‘‘.'ith the vrork of Rhyne et al. 

[22] v/ho used sputtered samples of TbFeg and showed that 
these samples are structurally amorpiaous and magnetic'll^ orde- 
red. They reported that botn the magxietization (li) and the 
critical temperature (1^,) were reduced in amorphous state. 

The reduction in magnetization c-nd To was attributed to the 
strong random local anisotropy [ 23 ]. later investigations 
[24-2 6] have established that the reduced Tc is d’ue to a 
reducedPe-Fe exchange interaction. Heiman and lee [27] from 
their studies on amorphous Lale found that in contrast to 
YFe and LuPe the sample were ferromagi'ietic with Tc necx room 
temperature. This suggested that the Fe-Fe exchange inter- 
action somehow depends on the size of the rare-earth because 
the most obvious difference betxxeen la, Y and Lu is orfly in 
the size. The radius of the R ahoms decreases witi, the atomic 
ntunber from 1.87 for la to 1.74 £ for lu (with b-ie exception 
of multiple- valence elements). The spin values and so the 
de C-ennes factor [IGF = (g-l)*^ J(J+l)] also varies systemati- 
cally, rising from la to Gd and then decreasing from Gd to Lu. 
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The systematic decrease of Tc with increasing atomic number 
is attributed to the decrease of the DIF. The extent to which 
each of the parameters, radius and DG-P, contribute to the 
systematic decrease of Tc has not yet been established. 

Recently Oullen [28] had introduced a model in which the 
magnetic interaction of the R atom aligns nearest neighbour 
Fe atoms in a spin-glass-like state. The work on amorphous 
YPe system by Coey et al. [29] suggests a local co-oi‘dination 
model in which the moment on the Fe atom is a function of 
the number of Fe nearest neighbours. 

Another interesting result obtained by Ghappert et al. 
[30] is the observation of transition metal moments larger in 
amorphous state than in crystalline counter part. A qualita- 
tive explanation for this observation is given in terms of the 
consequences of charge transfer from the rare-earth R to the 
empty d-orbitals of the transition metal [31]. The evidences 
available thus far suggests that all these models should be 
taken into account along with more experimental evidences for 
a comprehensive understanding of the amorphous HFe systems. 

Generally all these materials are ferromagnetic so that 
li is difference between the two sublattice magnetisations (R 
sublattice and Fe sublattice) and so they are liigiily composition 
dependent.. More over the spins with in a given satlottice 
are not necessarily aligned. In the case of Fe, the non 
alignment apparently results from large exchange fluctuations. 
The degree of non-alignment with in the two sublattices depend 
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upon the concentration and species of R, 

In many of the experiments reported above the amorphous 
R-Pe alloys were prepared by two-source co-evaporation of the 
elements and were deposited on substrate of fused-quartz or 
on polyparabanic acid sheets, A number of recent studies have 
demonstrated the potential of using an amorphous rapidly 
quenched alloy for the magnetic hardening of the rare-carth- 
iron or rare earth-iron based alloys [32-35] • In this vrork 
the magnetic and Mhssbauor studios on rapidly quenched amorphous 
^^0 33^^0 67 described. The properties of those alloys 

as well as those of the other light rarc-carths, are of particular 
interest because of their natural abundance makes them considera- 
bly less expensive, 

5 *3 *2 Experimental 

67 prepared by first molting in an 

induction furnace high purity Kd(99.9 ) and Fc( 99*99 y») 

(total 2 gm) undor high vacuum in a flat bortomed Zirconia 
crucible. The melting is repeated 2-3 times to get uniformity. 

Tho quenching was done by passing purified argon gas at low 
temperature through the molten alloy. The resulting brittle 
alloy was made into small pieces and annealed at 500 °0 for 
32 hours under high vacuum and quonchod in water. X-ray 
diffraction patterns of the powdered specimens were taken 
after taking precautions to prevent oxidations. Hie X-ray 
diffractrographs were characteristic of tho amorphous state 
(Pigs. 5.11 and 5.12). A rectangular piece of this material 
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Fig .512 X-ray diffraction pattern c* arrctphous Nd33Fe,^7 


0i loy 
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of size 0.214 X 0,294 X 0.152 cm was male out of this for 
magnetisation measurements in the PARC vibrating sample magne- 
tometer (VSM). Powdered sample of this alloy was mized with 
epoxy resin to form a thin plate for Mossbauer sfodies. 
M'dssbauer spectra were taken at different temperatures and 
the comparison with spectra taicen with normal powdered sample 
kept between celotape showed that the adhesive doe.-, not modify 
the spectrum. More over it prevents oxidation of the sample 
and ensures homogenity in distribution. 

Magnetisation measurements were carried out from room 
temperature to 425 K on the PARC vibrating sample magnetometer 
in fields up to 6.5 kOe. 

M’dssbauer effect measurements were made with the set up 
described in Chapter 2 in the constant acceleration mode from 
liquid nitrogen temperature to 475 *^1. 

5.3,3 Results and discussions 

(a) Results 

The temperature dependence of magnetisation of the 
amorphous Mq 5 ^ alloy is determined under residual 

field ( 40 gauss) during the increase and decrease of 

temperature. The results shows (Rig. 5.13) a sharp transition 
at around 333 'JRe thermol hysterisis present during 

heating and cooling is found to be negligibly small in this 
alloy. The Curie temperature (Tc) was obtained from the 
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da/dT vs T curve (Pig. 5.14) 8 .nd from the variation of with 

dT^ 

T. The value of Tc obtained in this way is 352,89 n. The 
temperature dependence of the inverse suscej)tibility is also 
shown in Pig. 5.13 and a linear dependence is observed. The 
Curie temperature of the recently synthesized EdPe'. compound 
[21] is higher than this value (Tc = 578 K) as erpected. 

Magnetisation measurements (o versus H) for 57 

alloy at several temperatures near Tc are also taiceu .:ith the 
VSM and are shown in Pig. 5.15. The temperature wrs hept stable 
within 0 . 2 *^, with the temperature controller developed by the 
author. It is observed that the sample is not getting saturated 
upto 6.5 hOe. The temperature dependence of magnetisation at 
various values of the applied field near the Curie temperature 
is shown in Pig. 5.16. The effective field acting on the 
sample is different from the applied field by a quantit;j^ called 
the demagnetising field of the sample. The demagneti sing field 
is given by = 4 mil . (magnetisation) where the magnetisation 
is in gauss and h is the geometrical demagnetisation fe.ctor 
which depends on the definite geometry of the sample. The 
geometric demagnetisation factor for a rectangular sanple was 
determined by using the expression [ 36 ] 

o . / 2,2, 2x1/2 

N = - Cot"^ [ia„±kJ-cJ_c^ 

•jt ■- a,D 

where the c axis is kept pareilel to the magnetic field. 
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The isotherms near the temperature Tc, drawn with 
H 

against — , called the Arrotts plots, are shown in iig, 5.17. 
Beyond 2 KGauss the Arrott plots are straight lines showing 
that the relation H/o = oc + p is very closely followed. The 
slope of these straight lines give ^ and the intercept on the 


a 

P 


a~ axis gives - ^ , where a and p are called the thernodynamic 


coefficients. The transition temperature is that where the 

thermodynamic coefficient a, changes sign and passes through 

zero (Big. 5. IS). The transition temperature (Tc) obtained 

by this method is 322 K, which is less than that obtained from 

the T curve by about 10 A. 

57 

The Fe Hdssbauer spectra at different temperatures of 


the amorphous 57 alloy is sho-ni in Fig, 5,19. The 

spectrum is typical of the amorphous structure where the line 
width goes on increasing as v/e go from the centre ox 'uhe 
spectrum. The spectra was recorded from 90 to 475 T* The 
data were analysed by a least square fit wizh Lorentzian line 
shape and also with the field distribution progrcnmies developed 
by the author. From 315 K up wards the spectrum shown the 
presence of both magnetic hyperfine splitting and qundrupole 
splitting superposed. The j.Tdssbauer parameters ob-cained at 
various temperature upto 320 K are shown in Table 5.5. The 
hyperfine field distribution P(H) clearly shows the presence 
of two distinct environments for the Fe ion with two different 
values for the mean field B. The average varies from 
317 ( 1 ) KOe at 90 K to 109(3) KOe at 320 F. and the corresponding 
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TABLE 5.5 

Observed Mbssbauer parajneters of the ferromagnetic phase 


of ^® 0.67 

amorphous 

alloy . 


Temperature 

K 

Heff 

Isomer Shift* 

(mm/ sec) 

Q* S . 

(mm/ sec ) 

90 

317(1) 

0.29 (1) 

0.27318) 

159 

279(2) 

0.074(9) 

0,356(9) 

168 

267 ( 2 ) 

0 . 152 ( 8 ) 

0.50 ( 1 ) 

216 

261(5) 

0.050 ( 1 ) 

0.26 ( 1 ) 

22 2 

251 ( 2 ) 

0 . 050 ( 1 ) 

0.24 (1) 

293 

173(3) 

-0.087(9) 

0.174(8) 

306 

155(2) 

-0.11 (1) 

0.26 (1) 

320 

109(3) 

- 0.13 ( 1 ) 

- 


*With respect to iron metal. 
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isomer shift change from +0,29(1) to -0.13(1 ) (with respect 
to iron metal) as the temperature is increased in the above 
range, figure 5.20 shows the variation of v/ith tempera- 
ture in NIq 67* value of ^^^^(O) obtained by 

extrapolation to zero degree Eelvin is 320 KOe and the value 
of To obtained from this curve for to be zero is 345 K. 

The hyperfine field distributions at various temperatures 
are shown in Pig, 5.21. The vaiues of H (mean fl.eld and the 
full width at half maxim-urn of the P(H) distribution ( *H) are 
shown in Table 5.6. 

In the paramagnetic phase amorphous alloy 

shows quadrupole splitting with well defined peal^s. liie 
quadrupole splitting at 333 E is 0.66 mm/sec with an isomer 
shift of -0.13 mm/sec with respect to j.ron metal. The quadru- 
pole splitting and the corresponding isomer shift obtained at 
various temperatures is shown in Table 5.7. The variation of 
quadrupole splitting with temperature obeys the empirical 
relation 

Q,S.(T) = K(1 - B T^/^) 

as can be seen from the Pig. 5.22. The value of the slope 

-5 .-“3/2 

pai’ameter obtained from tiiis graph is 6.93x10 K ' , 

The iron moment in the amorphous , 

^eff 

alloy using the relation by Gubbens et al. [37,35] is = I75~ 
is 2.17 Pg. 
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TABLE 5.6 

Temperature dependence of the most probable field 
mean fields (1) and EWM ( AH) of the hyperfine field 
distribution in M . Fe^ amorphous alloy. 


Temperature 

K 

^eff 

hOe 

*^0e 

ri2 

kOe 

kOe 

90 

317 

306 


45 

159 

279 

299 

2AG 

56 

168 

267 

299 

247 

42 

216 

261 

259 

225 

6C 

222 

251 

257 

2 23 

5C 

293 

174 

188 

14C 

35 

306 

155 

135 

- 

34 

315 

139 

128 

- 

48 
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TABLE 3.7 


Quadrupole splitting and isomer shift in the paramagnetic 
phase of Mq Ee^ amorphous alloy. 


Temperature 

K 

Q.S. 

(mm/sec) 


Isomer Shift* 
(mm/sec) 

525 

0.672 

(8) 


0.15 (1) 

333 

0,662 

(6) 


0.13 (1) 

343 

0.64 

(1) 


0.12 (1) 

363 

0.609 

(9) 


0.116(8) 

423 

0.53 

(1) 


0.11 (1) 

473 

0,45 

(1) 


- 


*¥ith respect to iron metal. 
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(b) Discussions 


The RPe 2 systems are generally ferrimagnets. The 
magnetisation measurements provide the global moment of a 
domain which is the sum of the net moments of the sublattices. 

On the other hand M'bssbauer measurements provides the net 
moment of the sublattice which depends on the magnitude and 
distribution in the direction of the individual atomic moments. 

In the crystalline phase H is aligned antiparallel to 
Pe, In the amorphous phase there is still antiparallel exchange 
between R and Pe, and there are also negative Pe-Pe bonds and the 
R-R coupling. Competion between these two .produces disorder 
and non-saturation and produces a random ferrimagnet or from 
Ooey's classification [29] a speromagnet (M’bssbauer spectra 
shows spatial isotropy in distribution). 

Another specific characteristic of tlae amorphous R~Pe 
alloys is that it shows a strong dependence on the size of R 
on Pig- 5.23 shows as a function of the atomic 

radius of the R-atom. The solid circles are data taken from 
Ref. [38]. The open circles are points for other R-Pe alloys. 

In these cases the contributions to is from the magnetic 
nature of the R- element such as core polarisation and dipolar 
terms because of the non zero value for the DG-P, Whereas for 
Lu, T and la the D&P = 0. Even then the trend is obeyed by 
these alloys also. 

The finite quadrupole interaction in the paramagnetic 
state of amorphous R-Pe systems are observed earlier also [ 39 ]* 
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Our measurements also shows the evidence for a non zero EPG- at 
the Pe sites. The loresence of a quadriipole interaction in the 
magnetic state shows strong local anisotropy at the Pe 'Site. 

If the material has a common ordering temperature for 
all ions and if the temperature dependence of the hyperfine 
field for all ions is identical then H should follow this 
temperature dependence and aH should decrease mo no tonic ally 
as one goes from low temperature to Tc. Tomiyoshi et al. [40] 
have shown that the thermal average <p> of the magnetic moments 
on which the weaker exchange field is acting decreases more 
rapidly with increase of temperature than the one with the 
stronger exchange fields. Therefore if the exchange field 
has local fluctuations at various lattice sites, the reduced 
magnetisation has a broad distribution at high temperature 
and also different ions may have hyperfine fields th:t show 
quite different temperature dependence and there will be a 
distribution of Curie temperatures in the syston which is 
observed in the case of the Mq 5 -y amorphous system. 

The distribution of EPG is evidenced by the p(lv!) 
distribution shown in Pig, 5.24. The distribution is broad 
near the transition temperature and goes on decreasing with 
increase of temperature since more and more Pe atoms comes 
to the paramagnetic state as the temperature is increased. 

5«3.4 Conclusion 

The magnetic properties of R— Pe alloys are quite varied 
and sensitive to the variation of temperature ahd other 
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parameters. M(3ssbauer spectroscopy can be used as a tool to 
probe into these amorphous systems to determine the individual 
magnetic moments. The study of an amorphous system also leads 
to a critical re-examination of the role of lattice in 
crystalline solids in the structural origin of hyperfine fields. 

It is hoped that the present work has contributed to 
our knowledge on the systematics of hyperfine interactions in 
metals, alloys and compounds and has helped bring us closer to 
the understanding of the origin of the hyperfine fields. 
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